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A review is presented of the development and current status of nuclear shell-model calculations 
in which the two-body effective interaction between the valence nucleons is derived from the free 
nucleon-nucleon potential. The significant progress made in this field within the last decade is 
Q ■ emphasized, in particular as regards the so-called T^ow-k approach to the renormalization of the 

^ . bare nucleon-nucleon interaction. In the last part of the review we first give a survey of realistic 

shell-model calculations from early to present days. Then, we report recent results for neutron-rich 
nuclei near doubly magic ^^^Sn and for the whole even-mass N = 82 isotonic chain. These illustrate 
^ . how shell-model effective interactions derived from modern nucleon-nucleon potentials are able to 

'sj" I provide an accurate description of nuclear structure properties. 
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1 Introduction 

The shell model is the basic framework for nuclear structure calculations in terms of nucleons. This 
model, which entered into nuclear physics more than fifty years ago [H E] , is based on the assumption 
that, as a first approximation, each nucleon inside the nucleus moves independently from the others in 
a spherically symmetric potential including a strong spin-orbit term. Within this approximation the 
nucleus is considered as an inert core, made up by shells filled up with neutrons and protons paired to 
angular momentum J = 0, plus a certain number of external nucleons, the "valence" nucleons. As is 
well known, this extreme single-particle shell model, supplemented by empirical coupling rules, proved 
very soon to be able to account for various nuclear properties [3], like the angular momentum and 
parity of the ground-states of odd-mass nuclei. It was clear from the beginning [4], however, that for 
a description of nuclei with two or more valence nucleons the "residual" two-body interaction between 
the valence nucleons had to be taken explicitly into account, the term residual meaning that part of the 
interaction which is not absorbed into the central potential. This removes the degeneracy of the states 
belonging to the same configuration and produces a mixing of different configurations. A fascinating 
account of the early stages of the nuclear shell model is given in the comprehensive review by Talmi [5]. 

In any shell-model calculation one has to start by defining a "model space", namely by specifying 
a set of active single-particle (SP) orbits. It is in this truncated Hilbert space that the Hamiltonian 
matrix has to be set up and diagonalized. A basic input, as mentioned above, is the residual interaction 
between valence nucleons. This is in reality a "model-space effective interaction", which differs from 
the interaction between free nucleons in various respects. In fact, besides being residual in the sense 
mentioned above, it must account for the configurations excluded from the model space. 

It goes without saying that a fundamental goal of nuclear physics is to understand the properties of 
nuclei starting from the forces between nucleons. Nowadays, the A nucleons in a nucleus are understood 
as non-relativistic particles interacting via a Hamiltonian consisting of two-body, three-body, and higher- 
body potentials, with the nucleon-nucleon (NN) term being the dominant one. 

In this context, it is worth emphasizing that there are two main lines of attack to attain this ambitious 
goal. The first one is comprised of the so-called ab initio calculations where nuclear properties, such as 
binding and excitation energies, are calculated directly from the first principles of quantum mechanics 
using appropriate computational scheme. To this category belong the Green's function Montecarlo 
Method (GFMC) P [7j, the no-core shell-model (NCSM)) [HE], and the coupled-cluster methods 
(CCM)) [TOl [11]. The GFMC calculations, on which we shall briefly comment in Sec. 12.31 are at 
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present limited to nuclei with A < 12. This limit may be overcome by using limited Hilbert spaces 
and introducing effective interactions, which is done in the NCSM and CCM. Actually, coupled-cluster 
calculations employing modern NN potentials have been recently performed for ^^O and its immediate 
neighbours [121 [131 [TH . 

The main feature of the NCSM is the use of a large, but finite number of harmonic oscillator 
basis states to diagonalize an effective Hamiltonian containing realistic two- and three- nucleon inter- 
actions [15j . In this way, no closed core is assumed and all nucleons in the nucleus are treated as active 
particles. Very recently, this approach has ben apphed to the study of A = 10 — 13 nuclei \16\. Clearly, 
all ab initio calculations need huge amount of computational resources and therefore, as mentioned 
above, are currently limited to light nuclei. 

The second line of attack is just the theme of this review paper, namely the use of the traditional 
shell-model with two-nucleon effective interactions derived from the bare A^A^ potential. In this case, 
only the valence nucleons are treated as active particles. However, as we shall discuss in detail in 
Sec. 131 core polarization effects are taken into account perturbatively in the derivation of the effective 
interaction. Of course, this approach allows to perform calculations for medium- and heavy-mass nuclei 
which are far beyond the reach of ab initio calculations. Many-body forces beyond the A^A^ potential 
may also play a role in the shell-model effective Hamiltonian. However, this is still an open problem 
and it is the main scope of this review to assess t he progress made in the derivation of a two-body 
effective interaction from the free A^A^ potential. 

While efforts in this direction started some forty years ago [T71 [TS] , for a long time there was 
widespread skepticism about the practical value of what had become known as "realistic shell model" 
calculations. This was mainly related to the highly complicated nature of the nucleon-nucleon force, in 
particular to the presence of a very strong repulsion at short distances, which made especially difficult 
solving the nuclear many-body problem. As a consequence, a major problem, and correspondingly a 
main source of uncertainty, in shell-model calculations has long been the two-body effective interaction 
between the valence nucleons. An early survey of the various approaches to this problem can be found 
in the Cargese lectures by Elliott |19], where a classification of the various categories of nuclear structure 
calculations was made by the number of free parameters in the two-body interaction being used in the 
calculation. 

Since the early 1950s through the mid 1990s hundreds of shell-model calculations have been carried 
out, many of them being very successful in describing a variety of nuclear structure phenomena. In the 
vast majority of these calculations either empirical effective interactions containing adjustable parame- 
ters have been used or the two body-matrix elements themselves have been treated as free parameters. 
The latter approach, which was pioneered by Talmi [20], is generally limited to relatively small spaces 
owing to the large number of parameters to be least-squares fitted to experimental data. 

Early calculations of this kind were performed for the p-shell nuclei by Cohen and Kurath [21], who 
determined fifteen matrix elements of the two-body interaction and two single-particle energies from 35 
experimental energies of nuclei from A = 8 up to A = 16. The results of these calculations turned out 
to be in very satisfactory agreement with experiment. 

A more recent and very successful application of this approach is that by Brown and Wildenthal 
[2^ [221 [21|, where the A = 17 — 39 nuclei were studied in the complete sd-shell space. In the final 
version of this study, which spanned about 10 years, 66 parameters were determined by a least squares 
fit to 447 binding and excitation energies. A measure of the quality of the results is given by the 
rms deviation, whose value turned out to be 185 keV. In this connection, it is worth mentioning that 
in the recent work of Brown and Richter [22] the determination of a new effective interaction for the 
s(i-shell has been pursued by the inclusion of an updated set of experimental data. As regards the 
empirical effective interactions, they may be schematically divided into two categories. The first is 
based on the use of simple potentials such as a Gaussian or Yukawa central force with various exchange 
operators consistent with those present in the interaction between free nucleons. These contain several 
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parameters, such as the strengths and ranges of the singlet and triplet interactions, which are usually 
determined by a fit to the spectroscopic properties under study. In the second category one may put 
the so-called "schematic interactions" . These contain few free parameters, typically one or two, at the 
price of being an oversimplified representation of the real potential. Into this category fall the well 
known pairing [261 1211 [28] , pairing plus quadrupole [29] , and surface delta [30l [311 [32] interactions. The 
successful spin and isospin dependent Migdal interaction [33] also belongs to this category. While these 
simple interactions are able to reproduce some specific nuclear properties [see, e.g., [34l[35]], they are 
clearly inadequate for detailed quantitative studies. 

We have given above a brief sketch of the various approaches to the determination of the shell-model 
effective interaction that have dominated the field for more than four decades. This we have done to 
place in its proper perspective the great progress made over the last decade by the more fundamental 
approach employing realistic effective interactions derived from modern A^A^ potentials. We refer to 
Ref. m for a detailed review of shell-model calculations based on the above empirical approaches. 

As we shall discuss in detail in the following sections, from the late 1970s on there has been sub- 
stantial progress toward a microscopic approach to nuclear structure calculations starting from the free 
A^A^ potential V^n. This has concerned both the two basic ingredients which come into play in this 
approach, namely the A^A^ potential and the many-body methods for deriving the model space effective 
interaction, Ves- As regards the first point, A^A^ potentials have been constructed, which reproduce 
quite accurately the A^A^ scattering data. As regards the derivation of Ves, the first problem one is 
confronted with is that all realistic A^A^ potentials have a strong repulsive core which prevents their 
direct use in nuclear structure calculations. As is well known, this difficulty can be overcome by resort- 
ing to the Brueckner G- matrix method (see Sec. 14.11) . which was originally designed for nuclear matter 
calculations. 

While in earlier calculations one had to make some approximations in calculating the G matrix for 
finite nuclei, the developments of improved techniques, as for instance the Tsai-Kuo method [Ml [37] 
(see Sec. 14.11) . allowed to calculate it in a practically exact way. Another major improvement consisted 
in the inclusion in the perturbative expansion for Ves of folded diagrams, whose important role had 
been recognized by many authors, within the framework of the so-called Q-box formulation ([38], see 
Sec.[SD. 

These improvements brought about a revival of interest in shell-model calculations with realistic 
effective interactions. This started in the early 1990s and continued to increase during the following 
years. Given the skepticism mentioned above, the main aim of this new generation of realistic calcula- 
tions was to give an answer to the key question of whether they could provide an accurate description 
of nuclear structure properties. By the end of the 1990s the answer to this question turned out to be in 
the affirmative. In fact, a substantial body of results for nuclei in various mass regions (see, for instance 
|39j ) proved the ability of shell-model calculations employing two-body matrix elements derived from 
modern A^A^ potentials to provide a description of nuclear structure properties at least as accurate 
as that provided by traditional, empirical interactions. It should be noted that in this approach the 
single- nucleon energies are generally taken from experiment (see, e.g.. Sec. 15.3.5]) . so that the calculation 
contains essentially no free parameters. Based on these results, in the last few years the use of realistic 
effective interactions has been rapidly gaining ground in nuclear structure theory. 

As will be discussed in detail in Sec. [SI the G matrix has been routinely used in practically all 
realistic shell- model calculations through 2000. However, the G matrix is model-space dependent as 
well as energy dependent; these dependences make its actual calculation rather involved (see Sec. 14. ip . 
In this connection, it may be recalled that an early criticism of the G-matrix method dates back to 
the 1960s [50]. which led to the development of a method for deriving directly from the phase shifts a 
set of matrix elements of V^n in oscillator wave functions [KJl [T^l E]. This resulted in the well-known 
Sussex interaction which has been used in several nuclear structure calculations. 

Recently, a new approach to overcome the difficulty posed by the strong short-range repulsion 
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contained in the free A^A^ potential has been proposed [l2l H3l HI]. The basic idea underlying this 
approach is inspired by the recent applications of effective field theory and renormalization group to 
low-energy nuclear systems. Starting from a realistic A^A^ potential, a low-momentum A^A^ potential, 
Viow-k, is constructed that preserves the physics of the original potential Vnn up to a certain cutoff 
momentum A. In particular, the scattering phase shifts and deuteron binding energy calculated from 
Vnn are reproduced by l^ow-k- This is achieved by integrating out, in the sense of the renormalization 
group, the high-momentum components of Vnn- The resulting Viow-k is a smooth potential that can 
be used directly in nuclear structure calculations without first calculating the G matrix. The practical 
value of the Viow~k approach has been assessed by several calculations, which have shown that it provides 
an advantageous alternative to the G-matrix one (see Sec I5.3.1"|) . 

The purpose of the present paper is to give a review of the basic formalism and the current status of 
realistic shell-model calculations and a self-contained survey of the major developments in the history 
of the field as regards both the A^A^ potential and the many-body approach to the derivation of the 
effective interaction. During the last four decades there have been several reviews focused on either of 
these two subjects and we shall have cause to refer to most of them in the following sections. Our review 
is similar in spirit to the one by Hjorth-jensen et al. |l5], in the sense that it aims at giving an overall 
view of the various aspects of realistic shell-model calculations, including recent selected results. The 
novelty of the present paper is that it covers the developments of the last decade which have brought 
these calculations into the mainstream of nuclear structure [see, for instance, |Sj]. As mentioned above, 
from the mid 1990s on there has been a growing success in explaining experimental data by means of 
two-body effective interactions derived from the free A^A^ potential, which has evidenced the practical 
value of realistic shell-model calculations. It is worth emphasizing that a major step in this direction has 
been the introduction of the low-momentum potential V^ow-k, which greatly simplifies the microscopic 
derivation of the shell-model effective interaction. On these grounds, we may consider that a first 
important phase in the microscopic approach to shell model, started more than 40 years ago, has been 
completed. It is just this consideration at the origin of the present review. 

Four more recent reviews [l6l [TTl WA HH] reporting on progress in shell-model studies are in some 
ways complementary to ours, in that they discuss aspects which we have considered to be beyond 
the scope of the present review. These regard, for instance, a phenomenologically oriented survey of 
shell-model applications [46J or large-scale shell- model calculations pT | Wl \ 

We start in Sec. [2] with a review of the A^A^ interaction, trying to give an idea of the long-standing, 
painstaking work that lies behind the development of the modern high-precision potentials. In Sec. [3] 
we discuss the derivation of the shell-model effective interaction within the framework of degenerate 
perturbation theory. The crucial role of folded diagrams is emphasized. Sec. H] is devoted to the 
handling of the short-range repulsion contained in the free A^A^ potential. We first discuss in Sec. 14.11 
the traditional Brueckner G-matrix method and then introduce in Sec. 14.21 the new approach based 
on the construction of a low-momentum A^A^ potential. In Sec. [5] we first give a survey of realistic 
shell-model calculations performed over the last four decades (Sees. 15.11 and . 15. 2p and then present 
some results of recent calculations. More precisely, in Sec. 15.3.11 a comparison is made between the 
G-matrix and \4ow-k approaches while in Sec. 15.3.21 results obtained with different A^A^ potentials are 
presented. In Sec. 15.3.31 we report selected results of calculations for nuclei neighboring doubly magic 
^^^Sn and compare them with experiment. Finally, in Sec. 15.3.41 we discuss the role of the many-body 
contributions to the effective interaction by investigating the results of a study of the even A^ = 82 
isotones The last section. Sec. [6l contains a brief summary and concluding remarks. 
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2 Nucleon-nucleon interaction 



2.1 Historical overview 

The nucleon-nucleon interaction has been extensively studied since the discovery of the neutron and in 
the course of time there have been a number of Conferences [l9l [501 EI| cind review papers [52l [53| 
marking the advances in the understanding of its nature. Here, we shall start by giving a brief historical 
account and a survey of the main aspects relevant to nuclear structure, the former serving the purpose 
to look back and recall how hard it has been making progress in this field. 

As is well known, the theory of nuclear forces started with the meson exchange idea introduced by 
Yukawa [55]. Following the discovery of the pion, in the 1950s many efforts were made to describe the 
nucleon-nucleon {NN) interaction in terms of pion-exchange models. However, while by the end of the 
1950s the one-pion exchange (OPE) had been experimentally established as the long-range part of Vnn, 
the calculations of the two-pion exchange were plagued by serious ambiguities. This led to several pion- 
theoretical potentials differing quite widely in the two-pion exchange effects. This unpleasant situation 
is well reflected in various review papers of the period of the 1950s, for instance the article by Phillips 
j56j : a comprehensive list of references can be found in Ref. [52] . 

While the theoretical efforts mentioned above were not very successful, a substantial progress in the 
experimental study of the properties of the A^A^ interaction was made during the course of the 1950s. 
In particular, from the examination of the pp scattering data at 340 MeV in the laboratory system 
[57j inferred the existence of a strong short-range repulsion, which he represented by a hard sphere for 
convenience in calculation. As we shall discuss in detail later, this feature, which prevents the direct 
use of Vnn in nuclear structure calculations, has been at the origin of the Brueckner G-matrix method 
(Sec. 14. ip and of the recent Mow-k approach (Sec. 14.21) . 

At this point it must be recalled that as early as 1941 an investigation of the possible types of 
nonrelativistic A^A^ interaction at most linear in the relative momentum p of the two nucleons and 
limited by invariance conditions was carried out by Eisenbud and Wigner [58]. It turned out that the 
general form of Vnn consists of central, spin-spin, tensor and spin-orbit terms. Some twenty years later, 
the most general Vnn when all powers of p are allowed was given by Okubo [59], which added a quadratic 
spin-orbit term. When sufficiently reliable phase-shift analyses of A^A^ scattering data became available 
(see for instance Ref. [BD]), these studies were a key guide for the construction of phenomenological 
A^A^ potentials. In the early stages of this approach, the inclusion of all the four types of interaction 
resulting from the study of Eisenbud and Wigner (1941), with the assumption of charge independence, 
led to the Gammel-Thaler potential [61] , which may be considered the first quantitative A^A^ potential. 
In this potential, following the suggestion of Jastrow (1951), a strong short-range repulsion represented 
by a hard core (infinite repulsion) at about 0.4 fm was used. As we shall see later, it took a decade 
before soft-core potentials were considered. 

In the early 1960s two vastly improved phenomenological potentials appeared, both going beyond 
the Eisenbud- Wigner form with addition of a quadratic spin-orbit term. These were developed by the 
Yale group [62] and by Hamada and Johnston (HJ) [63]. Both potentials have infinite repulsive cores 
and approach the one-pion-exchange-potential at large distances. Historically, the HJ potential occupies 
a special place in the field of microscopic nuclear structure. In fact, it was used in the mid 1960s in the 
work of Kuo and Brown [18j, which was the first successful attempt to derive the shell-model effective 
interaction from the free A^A^ potential. We therefore find it appropriate to summarize here its main 
features. This may also allow a comparison with the today's high-quality phenomenological potentials, 
as for instance Argonne Vig (see Sec. 12. 2p . The HJ potential has the form 

V = Vc{r) + VT{r)Su + VLs{r)L ■ S + VLL{r)Lu, (1) 
where C, T, LS and LL denote respectively central , tensor, spin-orbit and quadratic spin-orbit terms. 
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The operator 5*12 is the ordinary tensor operator and the quadratic spin-orbit operator is defined by 

L,, = [5Lj+{<T,-a,)]L^-{L.Sf. (2) 

The Vi (i=C, T, LS and LL) are spin-parity dependent, and hard cores, with a common radius of 
0.485 fm, are present in all states. With about 30 parameters the HJ potential model reproduced in a 
quantitative way the pp and np data below 315 MeV. 

As mentioned above, the era of soft-core potentials started in the late 1960s with the work of Reid 
[61] and Bressel et al. [65]. The original Reid soft-core potential Reid68 has been updated some 25 
years later [SB] producing a high-quality potential denoted as Reid93 (see Sec. 12. 2p . 

Let us now come back to the meson-theory based potentials. The discovery of heavy mesons in the 
early 1960s revived the field. This resulted in the development of various one-boson-exchange (OBE) 
potentials and in a renewed confidence in the theoretical approach to the study of the NN interaction. 
The optimistic view of the field brought about by the advances made during the 1960s is refiected in 
the Summary [19] of the 1967 International Conference on the Nucleon-Nucleon Interaction held at the 
University of Florida in Gainesville. A concise and clear account of the early OBE potentials (OBEP), 
including a list of relevant references, can be found in the review of the meson theory of nuclear forces 
by Machleidt [52]. 

During the 1960s sustained efforts were made to try to understand the properties of complex nuclei 
in terms of the fundamental A^A^ interaction. This brought in focus the problem of how to handle the 
serious difficulty resulting from the strong short-range repulsion contained in the free A^A^ potential. 
We shall discuss this point in detail in Sec. HI Here, it should be mentioned that the idea of overcoming 
the above difficulty by constructing a smooth, yet realistic, A^A^ potential that could be used directly 
in nuclear structure calculations was actively explored in the mid 1960s. This led to the development 
of a non-local, separable potential fitting two-nucleon scattering data with reasonable accuracy [HUES]- 
This potential, known as Tabakin potential, was used by the MIT group in several calculations of the 
structure of finite nuclei within the framework of the Hartree-Fock method [^ [701 ED 122] • An early 
account of the results of nuclear structure calculations using realistic A^A^ interactions was given at the 
above mentioned Gainesville Conference by Moszkowski [73]. 

As regards the experimental study of the A^A^ scattering, this was also actively pursued in the 
1960s (see [IH]), leading to the much improved phase-shift analysis of McGregor et al. [71], which 
included 2066 pp and np data up to 450 MeV. This set the stage for the theoretical efforts of the 1970s, 
which were addressed to the construction of a quantitative A^A^ potential (namely, able to reproduce 
with good accuracy all the known A^A^ scattering data) within the framework of the meson theory. 
In this context, a main goal was to go beyond the OBE model by taking into account multi- meson 
exchange, in particular the 27r-exchange contribution. These efforts were essentially based on two 
different approaches: dispersion relations and field theory. 

The work along these two lines, which went on for more than one decade, resulted eventually in the 
Paris potential [751 UHl HZl UHl [50] and in the so called "Bonn full model" [79], the latter including also 
contributions beyond 27r. In the sector of the OBE model a significant progress was made through the 
work of the Nijmegen group [HO]. This was based on Regge-pole theory and led to a quite sophisticated 
OBEP which is known as the Nijmegen78 potential. The Nijmegen, Paris, and Bonn potentials fitted the 
world A^A^ data below 300 MeV available in 1992 with a x^/datum = 5.12, 3.71, and 1.90, respectively 
[53]. 

To have a firsthand idea of the status of the theory of the A^A^ interaction around 1990 we refer 
to Ref. [SI] while a detailed discussion of the above three potentials can be found in [53]. Here we 
would like to emphasize that they mark the beginning of a new era in the field of nuclear forces and 
may be considered as the first generation of A^A^ realistic potentials. In particular, as will be discussed 
in Sec. 15.11 the Paris and Bonn potentials have played an important role in the revival of interest in 
nuclear structure calculations starting from the bare A^A^ interaction. We shall therefore give here a 
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brief outline of the main characteristics of these two potentials as well as of the energy-independent 
OBE parametrization of the Bonn full model, which has been generally employed in nuclear structure 
applications. 

In addition to the 27r-exchange contribution, the Paris potential contains the OPE and tu-meson ex- 
change. This gives the long-range and medium-range part of the NN interaction, while the short-range 
part is of purely phenomenological nature. In its final version [7S] the Paris potential is parametrized 
in an analytical form consisting of a regularized discrete superposition of Yukawa-type terms. This 
introduces a large number of free parameters, about 60 [53], that are determined by fitting the NN 
scattering data. 

As already mentioned, the Bonn full model is a field-theoretical meson-exchange model for the 
NN interaction. In addition to the 2'7r-exchange contribution, this model contains single vr, u, and 6 
exchanges and irp contributions. It has been shown [79j that the latter are essential for a quantitative 
description of the phase shifts in the lower partial waves while additional Sir and Air contributions are 
not very important. The Bonn full model has in all 12 parameters which are the coupling constants 
and cutoff masses of the meson-nucleon vertices involved. This model is an energy- dependent potential, 
which makes it inconvenient for application in nuclear structure calculations. Therefore, an energy 
independent one-boson parametrization of this potential has been developed within the framework of 
the relativistic three-dimensional Blanckenbecler-Sugar (BbS) reduction of the Bethe-Salpeter equation 
[79l[52]. This OBEP includes exchanges of two pseudoscalar (vr and rj), two scalar (cr and 6), and two 
vector (p and u) mesons. As in the Bonn full model, there are only twelve parameters which have to 
be determined through a fit of the NN scattering data. 

At this point, it must be pointed out that there are three variants of the above relativistic OBE 
potential, denoted by Bonn A, Bonn B and Bonn C. The parameters of these potentials and the 
predictions of Bonn B for the two-nucleon system are given in ^2] • The latter are very similar to the 
ones by the Bonn full model. The main difference between the three potentials is the strength of the 
tensor force as reflected in the predicted D-state probability of the deuteron P^. With Pd=4.4% Bonn 
A has the weakest tensor force. Bonn B and Bonn C predict 5% and 5.6%, respectively. Note that for 
the Paris potential Pd=5.8%. We shall have cause to come back to this important point later. 

We should now mention that there also exist three other variants of the OBE parametrization of the 
Bonn full model. These are formulated within the framework of the Thompson equation [52l [81] and 
uses the pseudovector coupling for vr and r], while the potential defined within the BbS equation uses 
the pseudoscalar coupling. It may be mentioned that the results obtained with the Thompson choice 
differ little from those obtained with the BbS reduction. A detailed discussion on this point is given in 

m 

As we shall see later, the potential with the weaker tensor force, namely Bonn A, has turned out 
to give the best results in nuclear structure calculations. Unless otherwise stated, in the following we 
shall denote by Bonn A, B, and C the three variants of the energy- independent approximation to the 
Bonn full model defined within the BbS equation. However, to avoid any confusion when consulting 
the literature on this subject, the reader may take a look at Tables A.l and A. 2 in [52] . 

2.2 High-precision NN potentials 

From the early 1990s on there has been much progress in the field of nuclear forces. In the first place, the 
A^A^ phase shift analysis was greatly improved by the Nijmegen group [82| |83| (HH [85]. They performed 
a multienergy partial-wave analysis of all A^A^ scattering data below 350 MeV laboratory energy after 
rejection of a rather large number of data (about 900 and 300 for the np and pp data, respectively) 
on the basis of statistical criteria. In this way, the final database consisted of 1787 pp and 2514 np 
data. The pp, np and combined pp + np analysis all yielded a x^/datum ^ 1, significantly lower than 
any previous multienergy partial-wave analysis. This analysis has paved the way to a new generation 
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of high-quality A^A^ potentials which, similar to the analysis, fit the A^A^ data with the almost perfect 
X^/datum ^ 1. These are the potentials constructed in the mid 1990s by the Nijmegen group, Nijml, 
Nijmll and Reid93 [SS], the Argonne Vig potential [HS], and the CD-Bonn potential [SZIEB]- 

The two potentials Nijml and Nijmll are based on the original Nijm78 potential (HO] discussed in 
the previous section. They are termed Reid-like potentials since, as is the case for the Reid68 potential 
|64j . each partial wave is parametrized independently. At very short distances these potentials are 
regularized by exponential form factors. The Reid93 potential is an updated version of the Reid68 
potential, where the singularities have been removed by including a dipole form factor. While the Nijmll 
and the Reid93 are totally local potentials, the Nijml contains momentum-dependent terms which in 
configuration space give rise to nonlocalities in the central force component. Except for the OPE tail, 
these potentials are purely phenomenological with a total of 41, 47 and 50 parameters for Nijml, Nijmll 
and Reid93, respectively. They all fit the A^A^ scattering data with an excellent x^/datum=1.03 [66] . 
As regards the D-state probability of the deuteron, this is practically the same for the three potentials, 
namely Pd in %= 5.66 for Nijml, 5.64 for Nijmll, and 5.70 for Reid93. It is worth mentioning that in 
the work by Stoks et al. p6j an improved version of the Nijm78 potential, dubbed Nijm93, was also 
presented, which with 15 parameters produced a x^/datum of 1.87. 

The CD-Bonn potential [88j is a charge-dependent OBE potential. It includes the vr, p, and uo mesons 
plus two effective scalar-isoscalar a bosons, the parameters of which are partial- wave dependent. As is 
the case for the early OBE Bonn potentials, CD-Bonn is a nonlocal potential. It predicts a deuteron 
D-state probability substantially lower than that yielded by the potentials of the Nijmegen family, 
namely Pd=4.85%. This may be traced to the nonlocalities contained in the tensor force [SB]. While 
the CD-Bonn potential reproduces important predictions by the Bonn full model, the additional fit 
freedom obtained by adjusting the parameters of the ai and a2 bosons in each partial wave produces 
a x^/datum of 1.02 for the 4301 data of the Nijmegen database, the total number of free parameters 
being 43. In this connection, it may be mentioned that the Nijmegen database has been updated [88] 
by adding the pp and np data between January 1993 and December 1999. This 1999 database contains 
2932 pp data and 3058 np data, namely 5990 data in total. The x^/datum for the CD-Bonn potential 
in regard to the latter database remains 1.02. 

The Argonne Vis model [86], so named for its operator content, is a purely phenomenological (except 
for the correct OPE tail) nonrelativistic A^A^ potential with a local operator structure. It is an updated 
version of the Argonne Vn potential [89], which was constructed in the early 1980s, with the addition 
of three charge-dependent and one charge-asymmetric operators. In operator form the Vig potential is 
written as a sum of 18 terms, 

= E (3) 

p=l,18 

To give an idea of the degree of sophistication reached by modern phenomenological potentials, it may 
be instructive to write here explicitly the operator structure of the Vi^ potential [86] . The first 14 charge 
independent operators are given by: 

Of=^'^^ = 1, (Tj ■ Tj), (cTj ■ CTj), (cTi • (Tj){Ti ■ Tj), Sij, Sij{Ti ■ Tj), L • S,L • S{Ti- Tj), 

L\ L\t, ■ T,), L\(T, ■ o-,), L\cTi ■ o-,)(t, ■ T,), {L . Sf, (L ■ 5)2(r, ■ r,). (4) 
The four additional operators breaking charge independence are given by 

where Tij = Sr^iTzj — ti ■ Tj is the isotensor operator analogous to the Sij operator. As is the case for 
the Nijml and Nijmll potentials, at very short distances the Vis potential is regularized by exponential 
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form factors. With 40 adjustable parameters this potential gives a x^/datum of 1.09 for the 4301 data 
of the Nijmegen database. As regards the deuteron D-state probability, this is Pd= 5.76%, very close 
to that predicted by the potentials of the Nijmegen family. 

All the high-precision A^A^ potentials described above have a large number of free parameters, say 
about 45, which is the price one has to pay to achieve a very accurate fit of the world A^A^ data. 
This makes it clear that, to date, high-quality potentials with an excellent x^/datum k, 1 can only be 
obtained within the framework of a substantially phenomenological approach. Since these potentials fit 
almost equally well the A^A^ data up to the inelastic threshold, their on-shell properties are essentially 
identical, namely they are phase-shift equivalent. In addition, they all predict almost identical deuteron 
observables (quadrupole moment and D / S-sidXe ratio) [51]. While they have also in common the 
inclusion of the OPE contribution, their off-shell behavior may be quite different. In fact, the short- 
range (high- momentum) components of these potentials are indeed quite different, as we shall discuss 
later in Sec. 14. 2[ This raises a central question of how much nuclear structure results may depend on 
the A^A^ potential one starts with. We shall consider this important point in Sec. 15.3.21 

The brief review of the A^A^ interaction given above has been mainly aimed at highlighting the 
progress made in this field over a period of about 50 years. As already pointed out in the Introduction, 
and as we shall discuss in detail in Sees. 15.11 and 15.21 this has been instrumental in paving the way to a 
more fundamental approach to nuclear structure calculations than the traditional, empirical one. It is 
clear, however, that from a first-principle point of view a substantial theoretical progress in the field of 
the A^A^ interaction is still in demand. It seems fair to say that this is not likely to be achieved along 
the lines of the traditional meson theory. Indeed, in the past few years efforts in this direction have 
been made within the framework of the chiral effective theory. The literature on this subject, which is 
still actively pursued, is by now very extensive and there are several comprehensive reviews ^SDl EB [S2] , 
to which we refer the reader. Therefore, in the next section we shall only give a brief survey focusing 
attention on chiral potentials which have been recently employed in nuclear structure calculations. 

2.3 Chiral potentials 

The approach to the A^A^ interaction based upon chiral effective field theory was started by Weinberg 
[93| [M] some fifteen years ago, and since then it has been developed by several authors. The basic 
idea [93] is to derive the A^A^ potential starting from the most general Lagrangian for low-energy pions 
and nucleons consistent with the symmetries of quantum chromodynamics (QCD), in particular the 
spontaneously broken chiral symmetry. All other particle types are "integrated out" , their effects being 
contained in the coefficients of the series of terms in the pion-nucleon Lagrangian. The chiral Lagrangian 
provides a perturbative framework for the derivation of the nucleon-nucleon potential. In fact, it was 
shown by Weinberg [9l] that a systematic expansion of the nuclear potential exists in powers of the small 
parameter Q/A^, where Q denotes a generic low-momentum and ^ 1 GeV is the chiral symmetry 
breaking scale. This perturbative low-energy theory is called chiral perturbation theory (xPT). The 
contribution of any diagram to the perturbation expansion is characterized by the power u of the 
momentum Q, and the expansion is organized by counting powers of Q. This procedure |9l] is referred 
to as power counting. 

Soon after the pioneering work by Weinberg, where only the lowest order A^A^ potential was obtained, 
Ordonez et a/. [95] extended the effective chiral potential to order (Q/A^)^ [next-to- next-to-leading order 
(NNLO), z^=3)] showing that this accounted, at least qualitatively, for the most relevant features of 
the nuclear potential. Later on, this approach was further pursued by Ordonez, Ray and van Kolck 
[nH [17], who derived at NNLO a A^A^ potential both in momentum and coordinate space. With 26 
free parameters this potential model gave a satisfactory description of the Nijmegen phase shifts up 
to about 100 MeV [97j. These initial achievements prompted extensive efforts to understand the A^A^ 
force within the framework of chiral effective field theory. 
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A clean test of chiral symmetry in the two-nucleon system was provided by the work of Kaiser, 
Brockmann and Weise [98j and Kaiser, Gerstendorfer and Weise [99j . Restricting themselves to the pe- 
ripheral nucleon-nucleon interaction, these authors obtained at NNLO, without adjustable parameters, 
an accurate description of the empirical phase shifts in the partial waves with L > 3 up to 350 MeV 
and up to about (50-80) MeV for the D-waves. 

Based on a modified Weinberg power counting, Epelbaum, Glockle and Meissner constructed a chiral 
NN potential at NNLO consisting of one- and two-pion exchange diagrams and contact interactions 
(which represent the short-range force) |100[ 1101] . The nine parameters related to the contact interac- 
tions were determined by a fit to the np S- and P-waves and the mixing parameter ei for -ELab < 100 
MeV. This potential gives a x^/datum for the A^A^ data of the 1999 database below 290 MeV laboratory 
energy of more than 20 [102] . 

In their program to develop a A^A^ potential based upon chiral effective theory, Entem and Machleidt 
set themselves the task to achieve an accuracy for the reproduction of the A^A^ data comparable to that 
of the high-precision potentials constructed in the 1990s, which have been discussed in Sec. 12.11 The 
first outcome of this program was a NNLO potential, called Idaho potential [103j . This model includes 
one- and two-pion exchange contributions up to chiral order three and contact terms up to order four. 
For the latter, partial wave dependent cutoff parameters are used, which introduces more parameters 
bringing the total number up to 46. This potential gives a x^/datum for the reproduction of the 1999 
np database up to ELab=210 MeV of 0.98 [TUl] . 

The next step taken by Entem and Machleidt was the investigation of the chiral 27r-exchange contri- 
butions to the A^A^ interaction at fourth order, which was based on the work by Kaiser |105l 1106] , who 
gave analytical results for these contributions in a form suitable for implementation in a next-to-next- 
to-next-to-leading (N^LO, fourth order) calculation. This eventually resulted in the first chiral A^A^ 
potential at N^LO [102] . This model includes 24 contact terms (24 parameters) which contribute to the 
partial waves with L < 2. With 29 parameters in all, it gives a x^/datum for the reproduction of the 
1999 np and pp data below 290 MeV of 1.10 and 1.50, respectively. The deuteron D-state probability 
is Pd= 4.51%. 

Very recently a A^A^ potential at N^LO has been constructed by Epelbaum, Glockle and Meissner 
(2005) which differs in various ways from that of Entem and Machleidt, as discussed in detail in [107] . 
It consists of one-, two- and three-pion exchanges and a set of 24 contact interactions. The total number 
of free parameters is 26. These have been determined by a combined fit to some nn and pp phase shifts 
from the Nijmegen analysis together with the nn scattering length. The description of the phase shifts 
and deuteron properties at N^LO turns out to be improved compared to that previously obtained by 
the same authors at NLO and NNLO [108] . As regards the deuteron D-state probability, this N^LO 
potential gives Pd=2.73-3.63%, a value which is significantly smaller than that predicted by any other 
modern A^A^ potential. 

In regard to potentials at N^LO, it is worth mentioning that it has been shown [ 109[ 1110] that the 
effects of three-pion exchange, which starts to contribute at this order, are very small and therefore of 
no practical relevance. Accordingly, they have been neglected in both the above studies. 

The foregoing discussion has all been focused on the two-nucleon force. The role of three-nucleon 
interactions in light nuclei has been, and is currently, actively investigated within the framework of ab 
initio approaches, such as the GFMC and the NCSM. Let us only remark here that in recent years the 
Green's function Monte Carlo method has proved to be a valuable tool for calculations of properties 
of light nuclei using realistic two-nucleon and three-nucleon potentials [6], llllj . In particular, the 
combination of the Argonne Vig potential and Illinois-2 three-nucleon potential has yielded good results 
for energies of nuclei up to ^^C |112] . For a review of the GFMC method and applications up to A=8 
we refer the interested reader to the paper by Pieper and Wiringa [7J. 

In this context, it should be pointed out that an important advantage of the chiral perturbation 
theory is that at NNLO and higher orders it generates three-nucleon forces. This has prompted appli- 
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cations of the complete chiral interaction at NNLO to the three- and four-nucleon systems |113j . These 
apphcations are currently being extended to light nuclei with A > 4 |114j . 

However, as regards the derivation of a realistic shell-model effective interaction the 3A^ forces have 
not been taken into account up to now. As mentioned in the Introduction, in this review we shall give 
a brief discussion of the three-body effects, as inferred from the study of many valence-nucleon systems 

3 Shell-model effective interaction 
3.1 Generalities 

As mentioned in the Introduction, a basic input to nuclear shell-model calculations is the model-space 
effective interaction. It is worth recalling that this interaction differs from the interaction between two 
free nucleons in several respects. In the first place, a large part of the A^A^ interaction is absorbed 
into the mean field which is due to the average interaction between the nucleons. In the second place, 
the A^A^ interaction in the nuclear medium is affected by the presence of the other nucleons; one has 
certainly to take into account the Pauli exclusion principle, which forbids two interacting nucleons to 
scatter into states occupied by other nucleons. Finally, the effective interaction has to account for 
effects of the configurations excluded from the model space. Ideally, the eigenvalues of the shell-model 
Hamiltonian in the model space should be a subset of the eigenvalues of the full nuclear Hamiltonian 
in the entire Hilbert space. 

In a microscopic approach this shell-model Hamiltonian may be constructed starting from a realistic 
A^A^ potential by means of many-body perturbation techniques. This approach has long been a central 
topic of nuclear theory. The following subsections are devoted to a detailed discussion of it. 

First, let us introduce the general formalism which is needed in the effective interaction theory. We 
would like to solve the Schrodinger equation for the A-nucleon system: 

=E,|^,), (6) 

where 

H = Ho + H,, (7) 

and 

A 

Ho = Y,i^^ + U{), (8) 

i<j=l i=l 

An auxiliary one-body potential Ui has been introduced in order to break up the nuclear Hamiltonian 
as the sum of a one-body term Hq, which describes the independent motion of the nucleons, and the 
interaction Hi. 

In the shell model, the nucleus is represented as an inert core plus n valence nucleons moving 
in a limited number of SP orbits above the closed core and interact through a model-space effective 
interaction. The valence or model space is defined in terms of the eigenvectors of Hq 

\<i>i) = [alal ... al]i\c), i = l,...,d, (10) 

where |c) represents the inert core and the subscripts 1, 2, n denote the SP valence states. The 
index i stands for all the other quantum numbers needed to specify the state. 
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The aim of the effective interaction theory is to reduce the eigenvalue problem of Eq. ([6]) to a 
model-space eigenvalue problem 

PH,sP\^^) = E^P\^^), (11) 

where the operator P, 

d 

P = 5^|$,)($,|, (12) 

i=l 

projects from the complete Hilbert space onto the model space. The operator Q = 1 — P is its 
complement. The projection operators P and Q satisfy the properties 

P^ = P, Q^ = Q, PQ = QP = . 

In the following, the concept of the effective interaction is introduced by a very general and simple 
method |115|, 1116] . Let us define the operators 

PHP = Hpp, PHQ = HpQ, 
QHP = Hqp, QHQ = Hqq . 
Then the Schrodinger equation ([6]) can be written as 

HppP\^,) + HpqQ\^,) = E,P\^,) (13) 
HqpP\^,)+HqqQ\^,) = E,Q\^,). (14) 



From the latter equation we obtain 



and substituting the r.h.s. of this equation into Eq. (IT^ we have 



= —^-—HqpP\^,), (15) 



Hpp + HpQ ^ Hqp ) P\^,) = E,P\^,) . (16) 



If the l.h.s. operator, which acts only within the model space, is denoted as 



Hea{Eu) = Hpp + HpQ— — — Hqp, (17) 

i^u - rlQQ 



Eq. (fT6!) reads 

i7eff(^.)P|^.) =P.P|^.), 



which is of the form of Eq. (JTTil . Moreover, since the operators P and Q commute with Hq, we can 
write Eq. (JT71) as 

H,s{E.) = PHoP + Ves{E,), (19) 

with 

V,siE,) = PH^P + PH,Q—^-—QH^P . (20) 

This equation defines the effective interaction as derived by Feshbach in nuclear reaction studies [116j . 
Now, on expanding {Ey — Hqq)~^ we can write 

^ ^ Ey — QHqQ ^ ^ Ey — QHqQ ^ Ey — QHqQ ^ ' ^ ^ 

13 




(22) 



Eqs. (120|) and ( l22l) are the desired result. In fact, they represent effective interactions which, used in 
a truncated model space, give a subset of the true eigenvalues. Bloch and Horowitz |115] have studied 
the analytic properties of the eigenvalue problem in terms of the effective interaction of Eq. (122|) . It 
should be noted, however, that the above effective interactions depend on the eigenvalue E^. This 
energy dependence is a serious drawback, since one has different Hamiltonians for different eigenvalues. 

Some forty years ago, the theoretical basis for an energy-independent effective Hamiltonian was set 
down by Brandow in the frame of a time- independent perturbative method [117] . Starting from the 
degenerate version of the Brillouin-Wigner perturbation theory, the energy terms were expanded out of 
the energy denominators. Then a rearrangement of the series was performed leading to a completely 
linked-cluster expansion. The energy dependence was eliminated by introducing a special type of 
diagrams, the so-called folded diagrams. 

A linked-cluster expansion for the shell- model effective interaction was also derived in Refs. [11811119] 
11201 1121] within the framework of the time-dependent perturbation theory. In the following subsection 
we shall describe in some detail the time-dependent perturbative approach by Kuo, Lee and Ratcliff 
|121] . We have tried to give a brief, self-contained presentation of this subject which as matter of 
fact is rather complex and multi-faceted. To this end, we have discussed the various elements entering 
this approach whithout going into the details of the proofs. Furthermore, we have found it useful to 
first introduce in Sees. 13.2.11 and 13.2.21 the concept of folded diagrams and the decomposition theorem, 
respectively, which are two basic tools for the derivation of the effective interaction, as is shown in 
Sec. 13.2.31 A complete review of this approach can be found in [38], to which we refer the reader for 
details. 

3.2 Degenerate time-dependent perturbation theory: folded-diagram ap- 
proach 

3.2.1 Folded diagrams 

In this section, we focus on the case of two- valence nucleons and therefore the nucleus is a doubly closed 
core plus two valence nucleons. We denote as active states those SP levels above the core which are 
made accessible to the two valence nucleons. The higher-energy SP levels and the filled ones in the core 
are called passive states. In such a frame the basis vector is 



In the complex time limit the time-development operator in the interaction representation is given 



$i) = K4]i|c) . 



(23) 



by 



[/(O, -cx))= lim lim U{0,t) = lim 



lim 

-oo(l— ie) 



,iHt-iHot 



e^0+ t^-oo(l-ie) e^0+ t- 



which can be expanded as 




y^HT / dt, / dt2... 

n=0 



dtr, H,{t,)H,{t2)...H,{Q, (24) 



where 
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Let us now act on with f/(0, — cxo): 

U{0, -oo)|$,) = f/(0, -oo)[a\al]i\c), 



(25) 



The action of the time-development operator on the unperturbed wave function may be represented 
by an infinite collection of diagrams. The type of diagrams we consider here is referred to as time-ordered 
Goldstone diagrams (for a description of Goldstone diagrams see for instance Ref. |117j ) 

As an example, we show in Fig. [1] a second-order time-ordered Goldstone diagram , which is one 
of the diagrams appearing in fl25|) . The dashed vertex lines denote VArAr-interactions (for the sake of 
simplicity we take Hi = V^n), a and b represent two passive particle states while 1, 2, 3, and 4 are 
valence states. From now on in the present section, the passive particle states will be represented by 
letters and dashed-dotted lines. 

a b 

I 1 

A A 



Fig. 1: A second-order time-ordered Goldstone diagram. 
The diagram A of Fig. [1] gives a contribution 



A = alal\c) X ^Kfe34'^^3412/(^), 

where Vajj^ys are non-antisymmetrized matrix elements of Vnn- I{A) is the time integral 



1 



(26) 



1(A) = lim lim (-if I dti 

e^0+ i'^-oo(l-ie) 



^j.^^-i(e:j+e4-ea — ei,)ti g— i(ei+e2-e3-e4)t2 



(27) 



where the e^'s are SP energies. A folded diagram arises upon factorization of diagram A, as shown 
in Fig. [21 From this figure, we see that diagram B represents a factorization of diagram A into the 
product of two independent diagrams. The time sequence for diagram A is > ti > t2 > t', while in 
diagram i? it is > ti > t' and > t2 > t', with no constraint on the relative ordering of ti and t2- 



1 f 2 
A 



t 1 



B 



1 t. 


4 a 


b 

A 






2 t, 







Fig. 2: A diagrammatic identity which defines the folded diagram. 
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Therefore, diagrams A and B are not equal, unless subtracting from B the time-incorrect contribution 
represented by the folded diagram C. It is worth noting that lines 3 and 4 in C are not hole lines, but 
folded active particle lines. From now on the folded lines will be denoted by drawing a little circle. 
Explicitly, diagram C is given by 



C 



where 



I{C) 



lim 

6^0+ t' 



lim (— i) 

— oo{l— ie) 



dti / dt2e 



-i{e3+e4- 



-<^b)tl g-i(ei+e2-e3-e4)t2 



(28) 



(29) 



Note that the rules to evaluate the folded diagrams are identical to those for standard Goldstone 
diagrams, except counting as hole lines the folded active lines in the energy denominator. 

We now introduce the concept of generalized folded diagram and derive a convenient method to 
compute it in a degenerate model space. Let us consider, for example, the diagrams shown in Fig. [31 
All the three diagrams A, B, and C have identical integrands and constant factors, the only difference 
being in the integration limits. The integration limits of diagram A correspond to the time ordering 
> ti > t2 > ts > t4 > t'. As pointed out before, the factorization of A into two independent diagrams 
(diagram B) violates the above time ordering. To correct for the time ordering in B, one has to subtract 
the folded diagram C, whose time constraints are > ti > t2 > t' , > > > t', and > ^2- Five 
different time sequences satisfy the above three constraints, thus C consists of five ordinary folded 
diagrams (see Fig. H]) and is called generalized folded diagram. From now on the integral sign will 
denote the generalized folding operation. 



j. 4 t 

t c d /> 

t 

3 4 



a b ' 



[ \h . . 

-) c d -ii X -ji a b -) 



3 4 

'i > '2 ' h > '4 



\ 

i 


i 

i 


\ c 


d i; 





A a b A 



'l > '2 . '3 > ^4 ' '3 > ' 2 



B 



Fig. 3: A diagrammatic identity which illustrates the generalized folded diagram. 

An advantageous method to evaluate generalized folded diagrams in a degenerate model space is 
as follows. Let us consider diagram C of Fig. [31 As pointed out before. A, B, and C have identical 
integrands and constant factors, so that we may write the time integral /(C) as ([38j, pp. 16-18) 



nC) = I{B) - I{A) 



ei + £2 - £3 - €4 [(es + €4 - Ce - e/)(e3 + £4 - - e^) 



(ei + 62 - ee - e/)(ei + £2 - - e^) 



ei + 62 - ea - efe 



(30) 



16 



In a degenerate model space, the first factor is infinite, while the second is zero. However, /(C) can 
be determined by a limiting procedure. If we write ei + €2 = €3 + €4 + A, where A — >• 0, we can put 
Eq. (I?T|) into the form 



lim — 

A^O A 



(es + €4- te- e/)(e3 + £4 - Cc - Cd) (es + £4 - Ce - e/ + A)(e3 + £4 - - 
1 



es + £4 - ea - ef, + A) 
1 



+ 



1 



(es + £4 - Ce - e/)2(e3 + 64 - Cc - td) (es + £4 - Ce - e/)(e3 + £4 - Cc - 
1 



(e3 + £4 - - eb) 

Thus, the energy denominator of the generalized folded diagram C can be expressed as the derivative of 
the energy denominator of the l.h. part of the diagram with respect to the energy variable u, calculated 
at ei + 62: 

m = ^ i '—^ . (32) 



£3 + e4 - - eb du \uj - - ed u - - tf J ^^^^^^^ 
The last equation will prove to be very useful to evaluate the folded diagrams. 




Fig. 4: Generalized folded diagram expressed as sum of ordinary folded diagrams. 



3.2.2 The decomposition theorem 

Let us consider again the wave function f/(0, — oo)|$i). We can rewrite it as 

t/(0, -oo)|$,) = f/L(0, -oo)|$,) X f/(0, -oo)|c), (33) 

where the subscript L indicates that all the Hi vertices in [7^(0, — oo)|$j) are valence linked, i.e. are 
linked directly or indirectly to at least one of the valence lines. We now factorize each of the two terms 
on the r.h.s. of Eq. (!33|) in order to write f/(0, — 00) in a form useful for the derivation of the 
model-space secular equation. 
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t = 

uq ( 0, -oo ) I c > ^ I c > + y y+ y y^. v + 

t =-oo 



(ClU (0,-00)10 = 1 + C)""""0 + 



+ 



+ 



Fig. 5: Diagrammatic representation of Uq{0, — oo)|c) and (c|t/(0, — oo)|c). 



We first consider f/(0, — cxd)|c), which can rewritten as 



f/(0, -oo)|c) = Uq{0, -oo)\c) X (c|f/(0, -oo)|c), 



(34) 



where Uq{0, — oo)|c) denotes the collection of diagrams in which every vertex is connected to the time 
t = boundary. In fact, Uq{0, —oo)\c) is proportional to the true ground-state wave function of the 
closed-shell system, while (c|?7(0, — oo)|c) represents all the vacuum fiuctuation diagrams. These two 
terms are illustrated in the first and second line of Fig. [5], respectively. 

A similar factorization of Ul{0, — oo)|<l>j) can be performed, which can be expressed in terms of the 
so-called Q-boxes. The Q-box, which should not be confused with the projection operator Q introduced 
in Sec. 13.11 is defined as the sum of all diagrams that have at least one ifi-vertex, are valence linked and 
irreducible (i.e., with at least one passive line between two successive vertices). Clearly, Ul{0, —oo)\^i) 
must terminate either in an active or passive state at t = 0, thus we can write 



[/L(0,-oo)|$.) = |xf) + |x?), 



(35) 



as shown in Fig. [61 Note that in this figure the intermediate indices k, a, ... represent summations over 
all P-space states. 



Ul(0,-oo) !(!>.) = £ 



j=i 



1 1 



J 

1 1 



+ I Ik + 



I I 
i 



n>d 



+ I la + 



1.1 
1 



Fig. 6: Diagrammatic representation of Ul{0, — oo)|$i). 

It is also possible to factorize out of \xf) a term belonging to \xf) by means of the folded-diagram 
factorization. In Fig. [7] we show, as an example, how a 2-Q-box sequence can be factorized. 

It is worth noting the similarity between Fig. [7] and Fig. [31 As a matter of fact, when factorizing 
diagram A time incorrect contributions arise in diagram B, that are compensated by subtracting from 
it the generalized folded diagram C. 

Using the generalized folding procedure, we are able to factorize out of each term in \xf) a diagram 
belonging to |xf ) (see Fig. [8]). Applying this factorization to all the terms in \xf), collecting columnwise 
the diagrams on the r.h.s. in Fig. [Hand adding them up, we may represent \xf) as shown in Fig. [91 
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t = n 




A 



Fig. 7: Factorization of a 2-Q-box. sequence. 

The collection of diagrams in the upper parenthesis of Fig. Ois simply ($j|f//^(0, — oo)|$i), which, 
according to Eq. (1551) . is related to |xf ) by 

d 

lxf) = El*^)(*^l^^(0'-^)l*^)- (36) 
Therefore, taking into account Figs. [6] and [9], and Eq. fl36|) we can express Ul{Q, — oo)|$j) as 

f/z.(0, -00)1$,) = Uql{0, -oo)|<l>,)(<l>,|t/z.(0, -oo)|<l>i), (37) 
where we have represented diagrammatically Uql{0, — oo)|$j) in Fig. [TOl 



n n j 




i 



Fig. 8: Factorization of \Xi) using the generalized folding procedure. 
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Eqs. fl33l) . fl34|) . and fl37j) define the decomposition theorem, which we rewrite as follows: 



where 
and 



f/(0, -oo)|<l>i) = J2 UQiO, -oo)|$,-)($,|f/(0, -oo)|$,), 



(<l>,-|f/(0,-oo)|$,) = (<l>,-|f/z.(0,-oo)|$,) X (c|f/(0,-oo)|c), 
[/q(0, -oo)|$,) = Uql{0, -oo)|$,) X Uq{0, -oo)|c) . 



(38) 

(39) 
(40) 
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Fig. 9: \xf) wave function. 

The decomposition theorem, as given by Eq. (1381) . states that the action of U{0, — oo) on can 
be represented as the sum of the wave functions Uq{0, —oo)\^j) weighted with the matrix elements 
($j|C/(0, — cx))|$j). Eq. (138!) will play a crucial role in the next Sec. I3.2.3[ where we shall give an 
expression for the shell-model effective Hamiltonian. 



J 

-co 





Fig. 10: Diagrammatic representation of Uql{0, —oc)\^j). 



3.2.3 The model-space secular equation 

For the sake of clarity, let us recall the model-space secular equation ( ITTi) 

where a = 1, ..,d, |\E'a) and are the true eigenvectors and eigenvalues of the full Hamiltonian H. 

From now on, we shall use, for the convenience of the proof, the Schrodinger representation. However, 
it is worth to point out that the results obtained hold equally well in the interaction picture. 
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First of all, we establish a one-to-one correspondence between some model-space parent states \pa) 
and d true eigenfunctions \^a)- Let us start with a trial parent state 

|pi) = ^El<^>.), (41) 



and act with the time development operator f/(0, — cxd) on it. More precisely, we construct the wave 
function 



(pi|f/(0, -oo)|pi) E^o+t'^-<x)(i~ie) (pi|e^^*'|pi) ■ 
By inserting a complete set of eigenstates of H between the time evolution operator and |pi), we obtain 



(pi|f/(0,-Oo)|pi) e^O+t' 

Here, is the lowest eigenstate of H for which (\E'i|pi) 7^ 0, this stems from the fact that the real 
exponential damping factor in the above equation suppresses all the other non-vanishing terms. 
This procedure can be easily continued, thus obtaining a set of wave functions 

I (pjt/(0,-oo)|p,)' ^ ^ 

where 

{pa\pa) = 1, (45) 

(Pal^a) 7^ 0, (46) 

(p„|^l) = (pj^2) = ... = (P«|^a-l) = . (47) 

The above correspondence (Eqs. H4lH7|) holds if the parent states \pa) are linearly independent. Under 
this assumption, we can write 



^.) = Y.ctm^ (4^ 



i=l 

with 



^ CfCf = . (49) 



i=l 



By construction, \^a) is an eigenstate of if, so, using Eq. (jH]), we can write 

^(0,-oo)|p^) ^ [/(0,-oo)|p^) 
(p„|[/(0,-oo)K) "(p„|[/(0,-oo)|p„) • ^ ^ 

Now, making use of Eq. ( 148|) and applying the decomposition theorem as expressed by Eq. ( l38l) . the 
above equation becomes 

E., ^q(0, -oo)|$,)($,|^(0, -oo)|$.)gr E., ^q(0. -oo)|$,)($,im -oo)l'^.)gf 

E.^(<f.|f/(o,-oo)|$™)C,"C;^ " Efc„.(<^'fc|f/(o,-oo)|<i>^)c,-c;^ • ^ ' 

In order to simplify the above expression, we define the coefficients 6" 

,a ^ Y..{^,\U{Q,-^)mCt ^ E.('^.I^L(0,-Oo)|$.)gr . . 

21 



where the r.h.s. of the above equation has been obtained by use of Eq. fl39l) . which cancels out the 
vacuum fluctuations diagrams {c\U{0, —oo)\c) of Fig. O Multiplying Eq. fIFIl) by it becomes 

d 

Y^i'^klHUgiO, -oo)|<l>,>J = E^bt, (53) 
j=i 

where use has been made of the relation ($fc|?7Q(0, —00)1$^) = 6kj (see Ref. |121] ). The above equation 
is the model-space secular equation we needed, where i^efr is given by HUq{0, — oo) and represents 
the projection of onto the model-space wave function 




Fig. 11: Goldstone linked-diagram expansion of Ec — E^. The cross insertions represent the Ho operator. 

In Eq. (1531) we can write the Hamiltonian H as H = Ho + Hi. First, let us consider the contribution 
from Ho- Since is an eigenstate of Ho and ($j|f/Q(0, — oo)|$j) = 6ij, we obtain 

mHoUgiO, -00)1$,) = {^i\Ho\^,) = 5,,{El + E°), (54) 

where E^ is the unperturbed core energy and E'^ is the unperturbed energy of the two valence nucleons 
with respect to E^. 

As for the matrix element (<l>j|iJit/Q(0, — oo)|$j), we see from inspection of Eq. (l40l) that it contains 
a collection of diagrams in which Hi is not linked to any valence line at t = 0. These diagrams 
are obtained acting with Hi on ^7^(0, — oo)|c) and their contribution to the l.h.s. of Eq. (1531) is 
5ij{c\HiUQ{Q, — oo)|c) = Sij{Ec — E^), Ec being the true ground-state energy of the closed-shell system. 
The diagram expansion of {E^. — E^) is given in Fig. [11] as illustrated in |122j . 

The other terms of ($j|iJif/Q(0, — oo)|$j) are all linked to the external active lines. By denoting, 
for simplicity, the collection of these terms as 

(<l>,|[//if/Q(0,-oo)]i|<l>,), (55) 

the secular equation (135]) can be rewritten in the following form: 

d 

K + Y.^H[HiUq{Q, -oo)]l|$,>J = {E^ - E,)bt . (56) 

We now define 



HIs = [HiUq{0,-oo)]l, (57) 

and show in Fig. [T2] a diagrammatic representation of its matrix elements, which has been obtained 
starting from the definition of Uql{0, — oo)|<l>j) given in Fig. [10] It should be noted that we have two 
kinds of Q-box, Q and Q'. Q and Q' are a collection of irreducible, valence-linked diagrams with at 
least one and two i/i-vertices, respectively. The fact that the lowest order term in Q' is of second order 
in Hi is just because of the presence of Hi in the matrix elements of Eq. flS^ . 
Formally, ifg^ can be written in operator form as 

Hls = Q-Q' Jq + Q' JqJq-Q' jQjQjQ+ ... , (58) 
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where the integral sign represents a generahzed folding operation. It is worth noting that, by definition, 
the Q-box contains diagrams at any order in Hi. Actually, when performing realistic shell- model 
calculations it is customary to include diagrams up to a finite order in Hi. A complete list of all the 
Q-hox diagrams up to third order can be found in Ref. [i5] . 



j j j k j 1 k 




i i k i 1 k i 



Fig. 12: Diagrammatic representation of H^g matrix elements. 

The Schrodinger equation is finally reduced to the model-space eigenvalue problem of Eq. 
whose eigenvalues are the energies of the A-nucleus relative to the core ground-state energy Ec- As 
mentioned above, the latter can be calculated by way of the Goldstone expansion [122j . expressed as 
the sum of diagrams shown in Fig. [TTl It is also worth noting that the operator H^^, as defined by 
Eq. ( l58|l . contains both one- and two-body contributions since in our derivation of the of Eq. ( l56l) we have 
considered nuclei with two- valence nucleons. All the one-body contributions, the so-called S-box |123j . 
once summed to the eigenvalues of Hq, E^, give the SP term of the effective shell-model Hamiltonian. 
The eigenvalues of this term represent the energies of the nucleus with one-valence nucleon relative to 
the core. This identification justifies the commonly used subtraction procedure [123] where only the 
two-body terms of H^^ (the effective two-body interaction V^g) are retained while the single-particle 
energies are taken from experiment. 

In this context, it is worth to mentioning that in some recent papers |124[ 1125] the SP energies and 
the two-body interaction employed in realistic shell-model calculations are derived consistently in the 
framework of the linked-cluster expansion. In particular, in Ref. |124j . where the light p-shell nuclei have 
been studied using the CD-Bonn potential renormalized through the V^ow-k procedure (see Sec. 14. 2p . 
a Hartree-Fock basis is derived which is used to calculate the binding energy of ^He and the effective 
shell-model Hamiltonian composed of one- and two-body terms. 

In concluding this brief discussion of Eq. (156|1 . it is worth to point out that for systems with more than 
two valence nucleons H^^ contains 1-, 2-, 3-,- ■ ■ , n-body components, even if the original Hamiltonian 
of Eq. ([7j) contains only a two-body force [126] . The role of these effective many-body forces as well as 
that of a genuine 3-body potential in the shell model is still an open problem (see [127] and references 
therein) and is outside the scope of this review. However, we shall come back to this point in Sec. 15.3.41 

In Sec. I3.2.1l we have shown how, in a degenerate model space, a generalized folding diagram can be 
evaluated. In [128] . it has been shown that a term like —Q' j Q may be written as 

-Q' Jq = ^^Q{u;), (59) 

where u is equal to the energy of the incoming particles. 

The above result may be extended to obtain a convenient prescription to calculate Hl^ as given by 
Eq. ( l58i) . We can write 

oo 
i=0 
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where 



Fo = Q, 
Fi = QiQ, 

F2 = Q2QQ + Q1Q1Q, (61) 



and 



1 d'^Qiu) 



(62) 



ujq being the energy of the incoming particles at t = — 00. 

Note that in Eq. (|6Tll we have made use of the fact that, by definition, 

dQ'{uo) dQ{uo) 



duj duj 



(63) 



The number of terms in grows dramatically with i. Two iteration methods to partially sum up 
the folded diagram series have been introduced in Refs. [371 1129i 1130] . These methods are known as the 
Krenciglowa-Kuo (KK) and the Lee-Suzuki (LS) procedure, respectively. In |130j . it has been shown 
that, when converging, the KK partial summation converges to those states with the largest model 
space overlap, while the LS one converges to the lowest states in energy. 

The LS iteration procedure was proposed within the framework of an approach to the construction of 
the effective interaction known as Lee-Suzuki method, which is based on the similarity transformation 
theory. In the following subsection we shall briefly present this method to illustate the LS iterative 
technique used to sum up the folded- diagram series (pUjl . 



3.3 The Lee-Suzuki method 

Let us start with the Schrodinger equation for the A-nucleon system as given in Eq. ([6]) and consider 
the similarity transformation 

H = X~^HX, (64) 

where X is a transformation operator defined in the whole Hilbert space. 
If we require that 

QHP = 0, (65) 

then it can be easily proved that the P-space effective Hamiltonian satisfying Eq. ffTTl) is just PTiP. 
Equation fl65|) is the so-called decoupling equation, whose solution leads to the determination of ifgff- 
There is, of course, more than one choice for the transformation operator X. We take 

X = e", (66) 

where the wave operator Q satifies the conditions: 

n = Qnp, (67) 

PnP = QVlQ = PVlQ = . (68) 
Taking into account Eq. fl67p . we can write X = 1 + Q and consequently 

H^g = PHP = PHP + PHiQn, (69) 
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with Hi defined in Eq. Qj, while the decouphng equation fl65|) becomes 

QHiP + QHQVL - VlPHP - nPHiQn = . (70) 

We now introduce the P-space effective interaction R by subtracting the unperturbed energy PHqP 
from Hes: 

R = PHiP + PHiQn . (71) 

In a degenerate model space PHqP = ujqP, we can consider a linearized iterative equation for the 
solution of the decoupling equation (!70!) 

[Eo - {QHQ - nr,-iPHiQ)]n^ = QHiP - Qn-iPH^P . (72) 

We now write the Q-box introduced in Sec. 13.2.21 in operatorial form as given in |130j 

Q = PHiP + PHiQ^—L--QHiP, (73) 

and define i?„ as the n-th order iterative effective interaction 



Rn = PHiP + PHiQVLr, . (74) 
fig = in Eq. fl72l) . i?„ can be written in terms of the Q"box and its 



Then, if we start with 
derivatives: 

Ri 

R2 

R3 



Q, 

[l-Qi]-'Q, 
[l-Qi-Q2R2]-'Q, 



-1 

Q (75) 



Rr. 



n-1 



n-1 



1 — Ql — ^ Qm Y\. 



m=2 fc=n— m+l 



where Qn is defined in Eq. fl62l) . 

The solution i?„ so obtained corresponds to a certain resummation of the folded diagrams to infinite 
order. In fact, if we consider, for example, R2 and expand it in power series of Qi, we obtain 

R2 = [l- QiY^Q = l + QiQ + QiQiQ + ... (76) 

It is clear from the above expression that R2 contains terms corresponding to an infinite number of 
folds. 

It can be shown that the Lee-Suzuki method yields converged results after a small number of 
iterations [123] , making this procedure very advantageous to sum up the folded-diagram series. 

4 Handling the short-range repulsion of the NN potential 

As already pointed out, a most important goal of nuclear shell-model theory is to derive the effective 
interaction between valence nucleons directly from the free NN potential. In Sec. El we have shown 
how this effective interaction may be calculated microscopically within the framework of a many-body 
theory. As is well known, however, Vnn is not suitable for this kind of approach. In fact, owing to the 
contribution from the repulsive core, the matrix elements of Vnn are generally very large and an order- 
by-order perturbative calculation of the effective interaction in terms of Vnn is clearly not meaningful. 
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Fig. 13: (Color online) Radial dependence of the relative wave function \E' of two nucleons interacting 
via a hard-core potential Vat at- ^ refers to the uncorrelated wave function. 



A resummation method has to be employed in order to take care of the strong short-range repulsion 
contained in V^n. 

This point is more evident if we consider the extreme model of a hard-core potential, as is the case 
of the potential models developed in the early 1960s. This situation is illustrated in Fig. [121 For these 
potentials the perturbation expansion of the effective interaction in terms of Vnn is meaningless, since 
each term of the series involving matrix elements of the NN potential between unperturbed two-body 
states is infinite. This is because the unperturbed wave function, in contrast to the true wave function, 
gives a non-zero probability of finding a particle located inside the hard-core distance. 

The traditional way out of this problem is the so-called Brueckner reaction matrix G, which is based 
on the idea of treating exactly the interaction between a given pair of nucleons |131j . The G matrix, 
defined as a sum of all ladder-type interactions (see Sec. 14.1.21) . is used to replace the NN interaction 
vertices once a rearrangement of the effective interaction perturbative series has been performed. 

Recently, a new method to renormalize the NN interaction has been proposed US]. A low- 
momentum model space defined up to a cutoff momentum A is introduced and an effective potential 
Viow-k is derived from Vnn- This Viow-k satisfies a decoupling condition between the low- and high- 
momentum spaces. Moreover, it is a smooth potential which preserves exactly the on-shell properties 
of the original potential and it is thus suitable to advantageously replace Vat a? in realistic many-body 
calculations. 

Sees. 14.11 and 14.21 are devoted to the description of the reaction G matrix and Vfow-k potential, 
respectively. 

Before doing this, however, it may be worth recalling that a method to avoid the G-matrix treatment 
to eliminate effects of the repulsive core in the A^A^ potential was proposed in the late 1960s [1321 BO] . 
As already mentioned in the Introduction, this consists in using the experimental A^A^ phase shifts to 
deduce matrix elements of the A^A^ potential in a basis of relative harmonic oscillator states. These 
matrix elements, which have become known as the Sussex matrix elements (SME), have been used 
in several nuclear structure calculations, but the agreement with experiment has been generally only 
semi-quantitative. A comparison between the results obtained by Sinatkas et al. |133] using the SME 
and those obtained with a realistic effective interaction derived from the Bonn A potential is made for 

[IMj. 



the N=50 isotones in 



4.1 The Brueckner G-matrix approach 
4.1.1 Historical introduction 



The concept of G matrix originates from the theory of multiple scattering of Watson |135l 1136] . In this 
approach, the elastic scattering of a fast particle by a nucleus was described by way of a transformed 
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potential obtained in terms of the Lippmann-Schwinger matrix [137] for the two-body scattering. The 
procedure of Watson for constructing such an "equivalent two-body potential" was generalized to the 
study of nuclear many-body systems by Brueckner and co-workers 1131^1138] . They introduced a reaction 
matrix for the scattering of two nucleons while they are moving in the nuclear medium. This matrix, 
which is known as the Brueckner reaction matrix, includes all two-particle correlations via summing all 
ladder-type interactions, and made it possible to perform Hartree-Fock self-consistent calculations for 
nuclear matter. 

Only a few years later Goldstone [122] proved a new perturbation formula for the ground-state energy 
of nuclear matter which gave the formal basis of the Brueckner theory. The Goldstone linked- diagram 
theory applies to systems with non-degenerate ground state, which is the case of nuclear matter as well 
as of closed-shell finite nuclei. 

The Brueckner theory was seen to be the key to solve the paradox on which the attention of many 
nuclear physicists was focused during the early 1950s. In fact, it allowed to reconcile a description of 
the nucleus in terms of an overall potential and the peculiar features of the two-body nuclear force. 
This was very well evidenced in the paper by Bethe [139] . whose main purpose was indeed to establish 
that the Brueckner theory provided the theoretical foundation for the shell model. 

After their first works, Brueckner and co-workers published a series of papers on the same subject 
[140[ I14H I142[ I143[ IT4:4j . where further analyses and developments of the method as well as numerical 
calculations were given. Important advances as regards nuclear matter were provided by the work of 
Bethe and Goldstone [145] and Bethe et al. [146] . For several years up to the 1960s, nuclear physicists 
were indeed very active in this field, as may be seen from the review papers by Day [147] , Rajaraman 
and Bethe [148] , and Baranger [149] , where comprehensive lists of references can be found. For a recent 
review of the developments in this area made Bethe and coworkers we refer to [150] . 

As regards open-shell nuclei, we have shown in the previous section that the model-space effective 
interaction may be obtained by way of a linked-diagram expansion containing both folded and non- 
folded Goldstone diagrams. An order-by-order perturbative calculation of such diagrams in terms of 
Vnn is not appropriate and one has to resort again to the reaction matrix G. 

The renormalization of Vnn through the Brueckner theory has been the standard procedure to derive 
realistic effective interactions since the pioneering works of the early 1960s. A survey of shell-model 
calculations employing the G matrix is given in Sees. 15.11 and 15.21 

In Sec. 14.1.21 we shall define the reaction matrix G and discuss the main problems related to its 
definition. Then, in Sec. 14.1.31 we shall address the problem of how the reaction matrix may be 
calculated, focusing attention on the Gt matrix for which plane waves are used as intermediate states. 
It is, in fact, this matrix which is most commonly used nowadays in realistic shell-model calculations. 
For simplicity, almost everywhere in this section we shall denote the NN potential by V. 

4.1.2 Essentials of the theory 

In the literature, the G matrix is typically introduced by way of the Goldstone expansion for the 
calculation of the ground-state energy in nuclear matter and closed shell nuclei. As a first step, the 
ground-state energy is written as a linked-cluster perturbation series. Then, all diagrams differing 
one from another only in the number of V interactions between two particles lines are summed. This 
corresponds to define a well-behaved two-body operator, the reaction matrix G, that replaces the 
potential V in the series. A very clear and simple presentation of the G matrix along this line is 
provided in the paper by Day [147] . 

Here, we shall not discuss the details of the G-matrix theory, but simply give the elements needed 
to make clear its definition in connection with the derivation of the model-space effective interaction 
for open-shell nuclei. Therefore, in the following we refer, as in the previous section, to an A-nucleon 
system with the Hamiltonian given by Eqs. ([7])- ([9]) and represented as a doubly closed core plus valence 
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nucleons moving in a limited number of SP orbits above the core. 
The two-body operator G is defined by the integral equation 



G{uj) = V + V- 



Q2p 



■G{uj) 



(77) 



where uj is an energy variable known as "starting energy" and an operator which projects onto 
particle-particle states, namely states composed of two SP levels above the doubly closed core. As we 
shall discuss later, this operator may be chosen in different ways depending on the specific context in 
which it is used. Here it is worth noting that its presence in Eq. ( 1771) reminds us that the G matrix, 
differently from the Lippmann-Schwinger T matrix, is defined in the nuclear medium. It may be also 
noted that the starting energy uj is not a free parameter. Rather, its value is determined by the physical 
problem being studied. For the T matrix, uo denotes the scattering energy for two particles in free space, 
but depends on the nuclear medium for the G matrix. In fact, we shall see that for a G matrix to be 
used within the linked-diagram expansion of the effective interaction uj depends on the diagram where 
G appears. Only in some cases it represents the energy of the two-particle incoming state, which was 
instead the meaning of the energy variable uj used in Sec. [31 
Let us start by writing the operator as 



Q2p = ^C{ij)\(f)ij){(f). 



(7J 



where the state \(f>ij), which is the antisymmetrized product of the two SP states \i) and is an 
eigenstate of Hq with energy ej + Cj. The constant C{ij) is 1 if the state pertains to the space 
defined by Q2p, otherwise. We now introduce the operator which projects on the complementary 
space 



2p 



1 ~ Q2p- 



Taking matrix elements of G between states of the P2p-space we have 

{(f)kl\G{uj)\(j)nm) = {(i)kl\V\(f)nm) + ^'^'^^ I ^1 '^^J ) 

which in a series expansion form becomes 

G{%2) 

'\(\)kl\G{uj)\<\)nm) = (0fcz|^|0nm) + y^(0H|^|0ii) {(pij\V\(pnm) 

'—^ UJ — ei — ej 



{(t)ij\G{u)\(i)nm) , 



(79) 



(80) 



OJ - ti 



U — tir — ej/ 



'\V\4>nm) + ■ 



The above expansion can be represented diagrammatically by the series shown in Fig. [TH The 
diagrams on the r.h.s. of this figure are known as "ladder" diagrams and each diagram corresponds to 
a situation in which a pair of particles interacts a certain number of times with the restriction that the 
intermediate states involved in the scattering must be those defined through the Q2p operator. In other 
words, two particles initially in a state of the P2p-space undergo a sequence of scatterings into states of 
the Q2p-space and then after several such scatterings go back to a state of the original space. 

For the sake of completeness, we introduce the correlated wave function lipnmi^^)) |145j . 



IV'nm(t^)) = |0nm) + 



UJ 



12) 
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Fig. 14: Diagrammatic representation defining the matrix G. 



which once iterated becomes 



IV'nm(w)) = \(f)nm) + 



Q2p 



uo-H, 



■G{^)\4'nm), 



^3) 



where use has been made of the integral equation (!77j) . Eq. (l83ll allows one to write the action of G on 
an unperturbed state as 

G'M|0„„) = F|^„„M), (84) 

which makes evident that the G operator may be considered an effective potential. The correlated wave 
function (182|) and its properties are extensively discussed in Refs. |149l 11511 1152j . 

At this point, we go further in our discussion clarifying the meaning of the starting energy u and 
illustrating some possible choices of the projection operator Q2p. 



Q2p=l 



Q2p=o 



n. n. 



Fig. 15: A graph of the projection operator Q2p appropriate for open-shell nuclei. 

Let us begin with Q2p. We define the Q2p operator by specifying its boundaries labelled by the three 
numbers (rii, 77,2,^3), each representing a SP level, the levels being numbered starting from the bottom 
of the potential well. Explicitly, using Eq. flTHl) the operator Q2p is written as 



1 ifi,j>n2 

1 ii rii < i < n2 and j > n2 

or vice- versa 
otherwise. 



The graph of Fig. fll5l) makes more clear its definition. Note that the index ns denotes the number 
of levels in the full space. In principle, it should be infinite; in practice, as we shall see in Sec. I4.1.3[ 
it is chosen to be a large but finite number. As regards the indices rii and ?t,2, the only mandatory 
requiriment is that none of them should be below the last occupied orbit of the doubly closed core. It 
is customary to take for ni the number the SP levels below the Fermi surface while n2 may be chosen 
starting from the last SP valence level and going up. How this choice is performed is better illustrated 
by an example. 

Let us consider the nucleus ^^O which in the shell-model framework is described as consisting of the 
doubly closed ^^O and two valence neutrons which are allowed to occupy the three levels of the sd shell. 
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The model-space effective interaction for ^^O may be derived by using the hnked-diagram expansion of 
Sec. [31 with the G matrix replacing the A^A^ potential in all the irreducible, valence- linked diagrams 
composing the Q-box. In so doing, one has to be careful to exclude from the Q-box those diagrams 
containing a ladder sequence already included in the G matrix. We may take the matrix G with a 
Q2p operator as specified, for instance, by (3, 6,00). In this case, one of two SP levels composing the 
intermediate state in the calculation of the G matrix has to be beyond the sd shell while the other one 
may be also an sd level. Then, when calculating the Q-box the included diagrams strictly depend on 
the considered G matrix. As an example, we have reported in Fig. [16] a first- and second-order diagram 
of the Q-box, the wavy lines denoting the G-interactions. As discussed in Sec. 13.2.21 the incoming and 
outcoming lines of the A and B diagrams are levels of the sd shell, while the the intermediate state of 
diagram B should have at least one passive line. This means that none of the two SP states, x or y, 
could be a Os or a Op level and at least one of them must be beyond the sd shell. Therefore, if we take 
the G matrix with the (3, 6, cxo) operator, all possible B diagrams of Fig. [THlare already contained 
in the A diagram. On the other hand, if we increase the n2 using a G matrix with the (3, 10, 00) 
Q2p operator, the B diagrams should explicitly appear in the calculation of the Q-hox providing that 
x,y <10 and x and/or y >7. A detailed discussion on the points illustrated here is in Ref. |153] . 



^\J'\J'\J'\J~\ X 



k Ik 1 

A B 



Fig. 16: Second-order ladder diagram contained in the Q-box. 

We now come to the energy variable uj which may be seen as the "starting energy" at which G is 
computed. We know that any G-vertex in a Q-box diagram is employed as the substitute for the ladder 
series of Fig. [TH and therefore the corresponding starting energy depends on the diagram as a whole 
and in particular on the location of the G-vertex in the diagram. For instance, the G matrices in the 
three vertices of diagrams A and B of Fig. [16] have to be calcutated at = + e/, corresponding to 
the energy of the two-particle incoming state. As an other example, we have shown in Fig. [T7| the 
diagram A. To illustrate how oj is evaluated we have written it as the sum of the ladder sequence Al, 
A2,...,An,.... The lowest vertex of diagram A corresponds to {4'mp\G{uj)\(l)kh) and its starting energy can 
be determined by looking at the two lowest vertices in diagram A2 whose contribution may be written 
as 

{<Pmp\V\(t)xy){(t)xy\V\<Pkh) ^^^-^ 
+ — ^x — ^y 

Note that the energy denominator of the G matrix is given hy u — Hq, Hq acting on the intermediate 
two-particle state \(pxy)- Therefore, to obtain the denominator of Eq. (!85|) as u — ex — e^, one has to 
choose u = ek + eh- Similarly, for the upper G- vertex of diagram A u = ek + ei + eh — em- In both the 
previous examples u is different from the energy of the incoming state, which results to be + e^. 

In concluding this subsection it is important to consider the dependence of the G matrix on the 
SP potential U, which explicitly appears in the denominator through Ho. The advantages of different 
choices have been extensively investigated essentially in connection with the Goldstone expansion and 
its convergence properties (see the review paper by Baranger [149j ). 
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Fig. 17: Illustration of the summation procedure to replace V vertices with G vertices. 



It is worth recalling that, as a general rule, in a perturbation approach to nuclear structure calcu- 
lations one starts by summing and subtracting the SP potential U to the Hamiltonian so as to have as 
perturbation Hi = V — U. As a consequence, the diagrams of the Q-hox for the shell-model effective 
interaction contain — f/ as well as V vertices. Diagrams where —U vertices are attached to particle lines 
between successive V interactions may be taken into account directly through the reaction matrix. This 
means that each diagram of the expansion shown in Fig. [H] is replaced with a new diagram in which 
any number of —U insertions has been introduced in the particle lines. A typical diagram is shown 
in Fig. [m In this way, a new reaction matrix may be introduced where the —U interactions for the 
intermediate particle states are summed up to all orders. More precisely, by changing the propagator 
of the particle lines from Hq to Hq — Q2pUQ2p one obtains the reaction matrix with plane waves as 
intermediate states, which is usually called Gt- Note that the term Q2pUQ2p instead of U is subtracted 
from Hq since, as pointed out in }149j , particle and hole states are defined with respect to the SP 
potential U . In other words, the Hamiltonian Hq — Q2pUQ2p, which commutes with Q2p, preserves the 
particle and hole states as defined by U. The reaction matrix Gt satisfies the integral equation 



Griuj) = V + VQ: 



2p 



UJ 



Q2ptQ2p 



Q2pGTiuj), 



(86) 



where one has to recall that t represents the kinetic energy. 

An extensive study of the properties and features of Gt in the framework of the shell model can be 
found in |153j . where arguments in favor of the use of Gt are given. We find it interesting to briefly 
recall these arguments here. 

First of all, the Gt matrix of Eq. (l86l) minimizes the dependence on the SP potential. In fact, it 
depends only on Q2p while the reaction matrix G of Eq. (177|) depends also on the spectrum of the SP 
states of the Q2p space. 
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X--- 



Fig. 18: Typical U insertions to the G- matrix intermediate states. 

The second argument is essentially based on the paper by Bethe et al. |146] . where it was shown 
that all —U and Brueckner-Hartree-Fock self-energy insertions cancel each other when the self-energy 
is calculated on energy shell. This result may be plausibly extended |153j to low-energy intermediate 
particle states not far off the energy shell, but such a cancellation is very unlikely for the high-energy 
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states. On the other hand, in nuclear matter calculations [152] it has been shown that the sum of 
self-energy-insertion diagrams and other three-body clusters is almost negligible, in other words they 
tend to cancel each other. If we assume that this nuclear-matter result still holds for finite nuclei, we 
reasonably expect that only the —U insertions survive for high-energy intermediate states, so leading 
to the reaction matrix Gt- The use of the G matrix of Eq. (177|) is instead suggested for low-lying 
states assuming that —U and Brueckner-Hartree-Fock self-energy insertions cancel approximately. In 
other words the reaction matrix calculation may be performed following a two-step procedure, which is 
discussed in detail in |153j . Here we only give a sketch of such procedure. 

Let us start by dividing the intermediate states of the Q2p space into high-energy and low-energy 
states, the Q2p operator being written as 

Q2p = Q2p + Q\p^ (87) 

with orthogonalized plane waves used for and harmonic oscillator wave functions for Q\.p. The fist 
step consists in the calculation of the Gt matrix of Eq. fl86l) with the projection operator taken as Q\^. 
Then using this reaction matrix instead of the potential the G matrix of Eq. fl77|) is calculated 
with a projection operator restricted to Q\^- Actually, the second step is carried out by means of a 
pertubative expansion which implies that the final reaction matrix, which we call Gt-, is written as 

G^{^) = Griuj) + GTiuj)-^hrGTiuj) + Gt{uj)-^^-Gt{uj)-^^-Gt{uj) + ■■■ (88) 

Hence Gt is obtained by summing to Gt the second -, third-, ■ ■ ■ order two-particle ladder diagrams with 
Gt{uj) vertices and intermediate states restricted to the low-lying states defined by the operator. 

These ladder terms may be directly included in the evaluation of the Q-box once the Gt of Eq. fl86l) with 
as projection operator has been calculated. Note that in so doing no —U or Brueckner-self-energy 

insertions are explicitly included in the Q-box diagrams since, on the basis of the previous arguments, 
one may assume that they cancel each other. 

It is worth to point out that no firm criteria exist on how to perform the partition of Eq. fl57|) . It 
is customary to take n2 in the definition of the Qgp operator so that it encompasses the shell above the 
model space |153] . For example, if we consider again ^^O we specify Q2p by (3, 10, oo) to define the 
Gt matrix, then while corresponds to the lpO/-shell and the ladder terms with two particle in this 
shell are taken into account directly when the Q-box is evaluated. 

4.1.3 Calculation of the reaction matrix 

Starting from the introduction of the G matrix, several approximate techniques have been developed 
for the solution of Eq. (177|) A review of them together with a critical discussion can be found in 
|147[ 11491 B5l 1154] . Among the early methods we would just like to mention the separation method of 
Moszkowski and Scott [15 5J and the reference spectrum method of Bethe et al. |146j . which have both 
been extensively applied to both nuclear matter and finite nuclei. For many years, they were the primary 
methods for G-matrix calculation and were first employed by Kuo and Brown fTB] to calculate from the 
Hamada- Johnston potential [63] the G matrix elements to be used in shell-model calculations for ^^O 
and ^^F. It is interesting to note that during the 1970s significant advances were made in the calculation 
of the reaction matrix as regards the G matrix of Eq. fl77|) with harmonic oscillator intermediate states 
as well as the Gt matrix of Eq. fISBl) having plane- wave intermediate states. In fact, accurate solutions 
of Eq. (177]) were proposed by Barrett et al. |156j and of Eq. fISB]) by Tsai and Kuo , both of them 
providing an exact treatment of the projection operator Q2p- 

Here we shall only focus on the technique to compute the reaction matrix of Eq. (!86|) . which we have 
used in our shell-model calculations. More precisely, we shall describe the method proposed by Tsai 
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and Kuo p6], which was further investigated and apphed by Krengiglowa et al. |153j . This method, 
based on a matrix inversion in the momentum space, allows, as mentioned above, an exact treatment of 
the projection operator except for the only approximation of using a finite value instead of infinity 
for ^3 defining Q2p. 

Let us start by noting that a solution of Eq. (186|) is 

Gt{u;) = V + VQ^p ^ , Q2pV, (89) 

Q2p[e - V)Q2p 

as it can be easily shown using the operator identity 

1111 

+ -jB—— , (90) 



A-B A A A-B 

with e = uj - t, A = Q2p e Q2p, and B = Q2pVQ2p- 

The propagator of Eq. ( l89l) was shown [36] to take the form 

^2p7; — 7 7jr77—Q2p = 77 7tP2p— — i — T^P'ip 77 , (91) 

Q2p{e - V)Q2p e-V e-V P2p-^P2p e-V 



where P2p is defined in Eq. f l79l) . The identity fl9T|) allows us to rewrite Eq. fl89|) as 

Gt{uj) = Gtf{uj) + AG{uj), (92) 

with 

^ -,._2_P,,_J_P,_4_V' (93) 

and 

Gtf(u;) = V + V^—V = V + V-Gtf{uj), (94) 
e — K e 

where we have used the identity fl90|) with = 1 to obtain relation fl94l) . 

An equivalent and very convenient expression of AG (a;) which has also the advantage to be applicable 
when V is a hard-core potential is 

AG{u) = -Gtf{u;)-P2p , ..p P2p-Gtf{uj). (95) 

e P2p[{l/e) + {l/e)GTF{^/e)\P2p e 

The set of Eqs. (!92l [94l [95l) provides a very convenient way to calculate the Gt matrix once a 
finite-ris approximation is made, namely a truncated P2p-space is taken. Let us define 

Gr{u) = GTF{l/e), 
GLiuj) = il/e)GTF, 
GLRiiu) = (1/e) + (l/e)GTF(l/e), (96) 

then Eq. (p2i) in a matrix formalism can be written as 

{(pkl\GT{uj)\(i)lm) = {(pkl\GTFiuj)\(f)lm) - ^ (0^ | Gl (t^) 1 0ij ) X {(f)ij\Gll^{uj)\(j)kl) {k l\G R{uj)\(f)lm) ■ (97) 

This equation makes clear the necessity of a finite- 713 approximation. In fact, the Gj^]^ matrix 
appearing in this equation is defined in the space and represents the inverse of the P2pGlrP2p 
matrix, which has to be a finite matrix. 
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The basic ingredient for the calculation of the Gt matrix through Eq. fl92l) is the matrix Gtf which 
is identical to the Gt defined by Eq. fl5^ . except that is free from the projection operator Q2p. A 
very convenient way to calculate Gtf{^^) is to employ the momentum-space matrix inversion method 
[157] ■ Once the Gtf matrix elements in the momentum space are known, the matrices Gl, Gr, and 
Glr can be easily evaluated in the same space. As a final step, these matrices are all transformed in 
the harmonic oscillator basis and the Gt matrix is obtained by means of Eq. (p7|l though simple matrix 
operations. 

Here we have only given a sketch of the procedure to solve Eqs. f l92l) . fl94|) . fl95l) . A detailed 
description of the whole procedure can be found in |153j . where the validity of an ria-finite approximation 
is discussed. It is worth mentioning that this point is also discussed in Ref. |45] . 

4.2 The Mow k approach 

In this subsection we shall describe in detail the derivation and main features of the low-momentum A^A^ 
interaction Mow-k, which has been introduced in Refs. |12[ HSl Hi] . The idea underlying the construction 
of Mow-k from a realistic model for Vnn is based on the renormalization group (RG) and the effective 
field theory (EFT) approaches |158l [92l [90} [9T] . As mentioned in Sec. [21 while the high-quality A^A^ 
potentials all reproduce the empirical deuteron properties and low-energy phase shifts very accurately, 
they differ significantly from one another in the high momentum behavior. A cutoff momentum A that 
separates fast and slow modes is then introduced and from the original Vnn an effective potential Viow-k, 
satisfying a decoupling condition between the low- and high-momentum spaces, is derived by integrating 
out the high-momentum components. The main result is that V^ow-k is free from the model-dependent 
high-momentum modes. In consequence, it is a smooth potential which preserves exactly the onshell 
properties of the original Vnn, and is suitable to be used directly in nuclear structure calculations. In 
Sec. 15.3.11 we will show that the Viow-k approach provides a real alternative to the Brueckner G-matrix 
renormalization procedure. 

4.2.1 Derivation of the low-momentum A^A^ potential 1 Viow-k 

In carrying out the above high-momentum integration, or decimation, an important requirement is that 
the low-energy physics of Vnn is exactly preserved by l^iow-k- For the two-nucleon problem, there is one 
bound state, namely the deuteron. Thus one must require that the deuteron properties given by Vat at 
are preserved by Viow-k- In the nuclear effective interaction theory, there are several well-developed 
model-space reduction methods. One of them is the Kuo-Lee-Ratcliff (KLR) folded diagram method 
[121[ [55]. which was originally formulated for discrete (bound state) problems. A detailed discussion of 
this method has been given earlier in Sec. [31 

For the two-nucleon problem, we want the effective interaction l^ow-k to preserve also the low-energy 
scattering phase shifts, in addition to the deuteron binding energy. Thus we need an effective interaction 
for scattering (unbound) problems. A convenient framework for this purpose is the T-matrix equivalence 
approach, as described below. 

We use a continuum model space specified by 



where p is the two-nucleon relative momentum and A the cutoff momentum which is also known as the 
decimation momentum. Its typical value is about 2 fm~^ as we shall discuss later. Our purpose is to 
look for an effective interaction PVcsP, with P defined above, which preserves certain properties of the 
full-space interaction Vnn for both bound and unbound states. This effective interaction is referred to 
as Viow-k- 




p <A 
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We start from the full-space half-on-shell T-matrix (written in a single partial wave channel) 



T{k\k,e) = VNN{k\k)+V / qHqVMN{k\q)- :^T{q,k,k' 

Jo - 



)■ 



This is the T-matrix for the two-nucleon problem with the Hamiltonian 

H = Hq + Vnn = t + Vnn , 
t being the relative kinetic energy. We then define a P-space low-momentum T-matrix by 



Tiow-k(p',P,P^) = Viow-k(p',p) + P / g^(ig\4ow-k(p',g)^ 5-Tiow-k 

Jo P - T 



(99) 



(100) 



(101) 



where {p',p) < A and the integration interval is from to A. In the above, the symbol V in front of 
the integration sign denotes principal value integration. We require the equivalence condition 



T{p',p,p^) = Tio^Mp',P,p'^); {p',p) < A. 



(102) 



The above equations define the effective low-momentum interaction; it is required to preserve the low- 
momentum (< A) half-on-shell T-matrix. Since phase shifts are given by the fuU-on-shell T-matrix 
T(p, p,p^), low-energy phase shifts given by the above Viow-k are clearly the same as those of V^n- 

In the following, let us show that a solution of the above equations may be found by way of the KLR 
folded-diagram method [12H [38] described in Sec. [3] as a means to construct the shell-model effective 
interaction. Within this approach the Viow-k may be written in the operator form as 



V, 



low— k 



Q-Q' Q + Q' Q Q-Q' Q Q Q + ... 



(103) 



which contains a Q-box whose explicit definition is given below. 

In the time-dependent formulation, the T-matrix of Eq. ( l99l) can be written as {k' \ VU (0, — oo) | k), 
U{0, — oo) being the time-evolution operator. In this way a diagrammatic analysis of the T-matrix can 
be made. A general term of T may be written as 



{k' 



V + v^—v + v^—v^—v + 



e(F 



e(P) e(F) 



k), 



(104) 



where e(fc^) = (A;^ - Hq). 

Note that the intermediate states (represented by 1 in the numerator) cover the entire space. In 
other words, we have 1 = P + Q where P denotes the model-space projection operator and Q its 
complement. Expanding it out in terms of P and Q, a typical term of T is of the form 



Q Q P Q P 
V-^V—V—V—V—V. 



(105) 



Let us now define the Q-box as 

{k' I Qi^) I k) = {k' 



Q Q Q 

V + v-^v + v-^v^^v + ... 



e(A;2 



e(P) e(F) 



A;), 



(106) 



where all intermediate states belong to Q. One readily sees that the P-space portion of the T-matrix 
can be regrouped in terms of a Q-box series, more explicitly as 



{P \ T\p) = {p 



. P . . P . P 

Q + Q-7-^Q +Q-7-^Q-7-^Q + 



e(p^ 



e(p^) e(p^) 



p). 



(107) 
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Fig. 19: Folded- diagram factorization of the half-on-shell T- matrix. 



Note that all the Q-boxes have the same energy variable, namely p^. This regrouping is shown 
in Fig. [T9|, where each Q-box is denoted by a circle and the solid line represents the propagator 
The diagrams A, B and C are respectively the one- and two- and three-Q-box terms of T, and clearly 
T=A-|-B-|-C-|-- ■ ■ . Note that the dashed vertical line is not a propagator; it is just a "ghost" line to 
indicate the external indices. 

Now a folded-diagram factorization for the T-matrix can be performed, following the same procedure 
as in Sec. [3] [ 12H [38] . As is discussed in Sec. 13.2. ![ diagram B of Fig. [19] may be factorized into the 
product of two parts (see Bl) where the time integrations of the two parts are independent from each 
other, each integrating from — oo to 0. In this way we have introduced a time-incorrect contribution 
which must be corrected. In other words, B is not equal to Bl, rather it is equal to Bl minus the 
folded-diagram correction B2. Schematically the above factorization is written as 

Q^Q = Q xQ-Q' J Q, (108) 

where the last term is the once-folded two Q-hox term. The meaning of the correction term B2 can also 
be visualized in terms of the differences in the Q-hox energy variables. With such energies explicitly 
written out, the above equation represents 

(p I Q(p')^)Q(p') I p) = Ilip' I W") I p") ip" I ^^(p') I - ^p' I ^' / ^ I p)- (109) 

Note that the energy variable for the first Q-box on the right hand side of this equation is p"^ instead 
of p^. Before factorization, the energy variable for this Q-box is p^. We recall that the integral sign 
represents a generalized folding operation and the leading Q-hox of any folded term must be at least of 
second order in Vjyjy. 

In the same way, we factorize the three-Q-box term C as shown in the third line of Fig. [THJ 

Q^Q^Q = Qx Q^Q-[Q' j Q] xQ + Q' J Q j Q, (110) 

Higher-order Q-box terms are also factorized following the same folded-diagram procedure. Let us 
now collect the terms in the figure in a "slanted" way. The sum of terms Al, B2, C3... is just the 
low-momentum effective interaction of Eq. (11031) . The sum Bl, C2, D3.... is Viow-k-Q. Similarly the 
sum Cl-t-D2-t-E3+- • • is just Viow-k—Q—Q- It is worth to mention that diagrams Dl, D2, ■ ■ ■ , El, E2, 
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• • • are not shown in Fig. [191 Continuing this procedure, it can be seen that the P-space portion of 
T-matrix is ^ow-k + ^low-k-f^". Namely, the Viow-k defined by Eq. fll03p is a solution of Eqs. fl^ - fll02p . 

An advantageous method to calculate the series in Eq. fllOSp is the Lee-Suzuki one, as described in 
Sec. 13.31 Following this algebraic approach the l^ow-k is expressed as 

liow-k = PVnnP + PVnnQuj, (111) 

where the operator u is given by the solution of the decoupling equation fl70p . Note that in this latter 
equation the letter Q was used instead of u to avoid confusion with the energy variable u defined in the 
same section. Here, since there is no problem of ambiguity, the same letter of the original paper by Lee 
and Suzuki (129j has been used and, for the sake of clarity, Eq. f lTOj) is re-written for the Hamiltonian 

(mil 

QVnnP + QHQuJ -ujPHP -ujPVnnQ^ = ^ ■ (112) 

The iterative technique presented in Sec. 13.31 to solve the decoupling equation is formulated for 
degenerate model spaces. Therefore, it is not useful in this case since PHqP = PtP is obviously non- 
degenerate. We can then resort to the iterative technique for non-degenerate model spaces proposed 
by Andreozzi |159j to solve the decoupling equation [see Eq. ( 11120 ] of the Lee-Suzuki method. This 
procedure, which will be referred to the Andreozzi-Lee-Suzuki (ALS) procedure, is quite convenient for 
obtaining the low-momentum NN interaction. 

We give now a sketch of this procedure. Let us define the operators: 



p{u) = PHP + PHQuj, (113) 
q{uj) = QHQ-uPHQ, (114) 



in terms of which the basic equations of the ALS iterative procedure read 



Xo = -{QHQ)-'QHP, 
xi = g(xo)~^xop(xo), 

Xn = g(xo + xi + ... + x„__i)"-^x„_ip(xo + Xi + ... + x„_i). (115) 
Once the iterative procedure has converged, x„ — >• 0, [159] . the operator u is given by 

n 

UJn = '^Xi, (116) 
i=Q 

and the V^ow-k is obtained by Eq. fillip . 

In applying the ALS method, we have employed a momentum-space discretization procedure by 
introducing an adequate number of Gaussian mesh points |153] . 

4.2.2 Phase-shift equivalent Hermitian V^iow-k's 

It should be pointed out that the low-momentum A^A^ interaction given by the above folded-diagram 
expansion is not Hermitian. This is not convenient for applications, and one would like to have an 
interaction which is Hermitian. We shall describe how a family of Hermitian l^ow-k's can be derived 
if we relax the half-on-shell T-matrix preservation. This subject has been studied in jl60] . where a 
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general framework for constructing hermitian effective interactions is introduced, and it is shown how 
this generates a family of hermitian phase-shift equivalent low-momentum A^A^ interactions. 

In the following, we outline the essential steps for the derivation of a Hermitian Vfow-k- We denote 
the folded-diagram low- momentum NN interaction obtained from the ALS procedure as Vls- The 
model-space secular equation for Vls is 

P{Ho + VLs)P\Xm) = ^mlXm), (117) 

where {Em} is a subset of the eigenvalues {En} of the full-space Schrodinger equation, {Hq + = 
Enl'^n), and \xm) is the P-space projection of the full-space wave function |\l/m), namely \xm) = -P|^m)- 
The above effective interaction may be rewritten in terms of the wave operator uj [see Eqs.( l66i) .(l671). 



PVlsP = Pe-'^iHo + V)e'^P - PHoP. (118) 

While the full-space eigenvectors |\1'„) are orthogonal to each other, the model-space eigenvectors 
|Xm)'s are clearly not so and consequently the effective interaction Vls is not Hermitian. We now make 
a Z transformation such that 

{Vm\Vm') = Smm'] m,m' = l,d, (119) 

where d is the dimension of the model space. This transformation reorients the vectors |Xm)'s such 
that they become orthonormal. We assume that the |Xm)'s are linearly independent so that exists, 
otherwise the above transformation is not possible. Since \vm) and Z exist entirely within the model 
space, we can write \vm) = P\vm) and Z = PZP. 
Using Eq. flllQp . we transform Eq. (11171) into 

Z{Ho + VLs)Z-'\vm) = Em\vm), (120) 

which implies 

Z{Ho + Vls)Z-^ = I Vm){Vm \ • (121) 

meP 

Since E^ is real (it is an eigenvalue of (Hq+V) which is Hermitian) and the vectors \vm) are orthonormal, 
Z{Hq + Vls)Z~^ is clearly Hermitian. The non-Hermitian secular equation (11171) is now transformed 
into a Hermitian model-space eigenvalue problem 

P{Ho + Vherm)P\Vm) = Em\Vm) , (122) 

with the Hermitian effective interaction given by 

Vherm = Z{Ho + Vls)Z-' - PH^P, (123) 

or equivalently 

Vherm = Zc-'^^H^ + V)e'' Z'^ - PHqP (124) 

To calculate Vherm, we must first have the Z transformation. Since there are certainly many ways 
to construct Z, this generates a family of Hermitian effective interactions, all originating from Vls- For 
example, we can construct Z using the familiar Schmidt orthogonalization procedure. We denote the 
hermitian effective interaction obtained using this Z transformation as Kchm- 

We can also construct Z analytically in terms of the wave operator. From the properties of the wave 
operator u, we have 

{Xm\ {I +UJ^Uj) \ Xm') = Smm'- (125) 
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It follows that an analytic choice for the Z transformation is 

Z = P(l + cu+Lu)i/2p_ ^126) 

This leads to the Hermitian effective interaction 

Kkb = P(l + cc.+u;)'/'P(i/o + + uj^uj)-^l''P - PH^P, (127) 

or equivalently 

Kkb = P{1 + uj^uj)-^'^{l + uj+){Ho + V){1 + uj){l + uj+uo)-^'^P - PHoP. (128) 

This is just the Okubo Hermitian effective interaction [161] . The Hermitian effective interaction of 
Suzuki and Okamoto |162l [TU] is equivalent to the above. 

There is another interesting choice for the transformation Z. As pointed out in |159j . the positive 
definite operator P(l + u!^uj)P admits the Cholesky decomposition, namely 

P{l + uj+uj)P = PLL'^P, (129) 

where L is a lower triangular matrix, being its transpose. Since L is real and within the P-space, 
we have 

Z = L^, (130) 
and the corresponding Hermitian effective interaction from Eq. f ll23p is 

Khoi = PL^P{Ho + Vls)P{L-YP PHoP (131) 

We have found that this Hermitization procedure is numerically more convenient than the other meth- 
ods. In the practical applications presented in this review we have employed this procedure. From now 
on, we shall denote the Hermitian low-momentum A^A^ interaction as V|ow-k- 

Using a solvable matrix model, it has been shown that the above Hermitian effective interactions 
V^chm, Vokh, and V^hoi can be in general quite different |160j from each other and from Vls, especially when 
Vls is largely non-Hermitian. For the Vat at case, it is fortunate that the Viow-k corresponding to Vls 
is only slightly non-Hermitian. As a result, the Hermitian low-momentum A^A^ interactions V^chm; Vokh, 
and V^hoi are all quite similar to each other and to the one corresponding to Vls, as discussed in |160] . 

We have just shown that a family of Hermitian effective interactions Viow-k's can be generated 
starting from the Lee-Suzuki transformation. In the following, we shall prove that all such interactions 
preserve the fuU-on-shell T- matrix T{p, p, p^), p < A, and consequently they are all phase-shift equivalent 

Let us consider two T-matrices 



T,{p') = V, + VMp')Tiip), (132) 

and 

T2{p'') = V2 + V2g{p'')T2{p'), (133) 

which are defined starting from two different potentials, Vi and V2, and/or two different propagators, 
(yf(p^) = V{pi — Hq)~^ and glp'"^) = V{p'2 — Hq)~^. These T-matrices are related by the well known 
two-potential formula [146] 

tKp") = Up') + nip"Mp") - 9ip')mp') + X2{p")[V2 - V^]X^ip'), (134) 
where the operator X = e'^ [see Eq. (!66ll ] is defined so as to have T = VX. 
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Table 1: Calculated binding energies BE^ of the deuteron (MeV) for different values of the cutoff A 
(fm~^). The V^ow-k has been calculated from the CD-Bonn potential for which BE^ = 2.224575 MeV 



A 0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


2.2246 


2.2246 


2.2246 


2.2246 


2.2246 


2.2246 



Applying the above relation to the half-on-shell T-matrix elements in momentum space, we have 

{p'\tKp")\p) = {p'\Tiip')\p) + {p'\tKp"Mp") - 9ip')]Tiip')\p) + {Mp')m - Vi)\Mp))- (135) 

Here the true and unperturbed wave functions are related by \ipi{p)) = Xi{p'^)\p) and similarly for ip2- 
Using the above relation, we shall now show that the phase shifts of the full-space interaction V^n are 
preserved by the Hermitian interaction Viow-k for momentum < A. Let us denote the last term of Eq. 
(11351) as D{p',p). We use Viow-k for V2 and Vnn for Vi. Recall that the eigenfunction of {Hq + Viow-k) 
is \vm) [see Eq. (11171) ] and that for H = {Hq + V^n) is l^&m). We define a wave operator 

Up=J2M{^rn\. (136) 

Then \vm) = f^pl^&m) and P\4ow-k-P = Up{Hq + Vnn)Up — PHqP. For our present case, {ip2{p')\ is {vp'\ 
and \ipi{p)) is l^&p). Since {vp'\Up = we have 

Dip',p) = i^p^KHUj, - Uj,H)\^iIp) = ip"~p')M^p). (137) 

Clearly D{p,p) = 0. The second term on the right hand side of Eq. (11351) vanishes when p' = p. Hence 

(p|Tiow-k(p')b) = {p\T{p^)\p), P<A, (138) 

where Tiow-k is the T-matrix for {PHqP + Viow-k) and T for {Hq + V^jy). Thus we have shown that 
the phase shifts of V^n are preserved by Viow-k- Recall that our T-matrices are real, because of the 
principal-value boundary conditions employed. 

4.2.3 Properties of Mow-k 

We have shown above that low-energy observables given by the full-momentum potential V^n are 
preserved by the low-momentum potential V^ow-k- This preservation is an important point, and it is 
worth checking numerically that the deuteron binding energy and low-energy phase shifts given by VivAr 
are indeed reproduced by V^ow-k- In this subsection we would like to report some results in this regard, 
together with some discussion about the differences between the various modern potentials. 

Before presenting the results, let us first address an important question concerning the choice of the 
momentum cutoff A in the derivation of Viow-k- Phase shifts are given by the full-on-shell T-matrix 
T(p, p,p^), hence, for a chosen A, Viow-k can reproduce phase shifts up to Eiab = 2ti?h? /M, M being the 
nucleon mass. Realistic NN potentials are constructed to fit empirical phase shifts up to Eiab ~ 350 
MeV [85j, which is the pion production threshold. It is reasonable then to require Viow-k to reproduce 
phase shifts also up to this energy. Thus one should use A in the vicinity of 2 fm~^. 

We first check the deuteron binding energy BE^ given by Viow-k- For a range of A, such as 0.5 
fm~^ < A < 3 fm~^, BEd given by V^ow-k agrees very accurately (to 4 places after the decimal point) 
with that given by Vnn- In Tabled! we present the results for the Viow-k derived from the CD-Bonn 
NN potential [88j. 
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We have also checked the phase shifts and the half-on-shell T{p\p,p'^) matrix elements with {p',p) < 

A. 

In Table [2] we report the neutron-proton ^5*0 phase shifts calculated both with the CD-Bonn potential 
and its l^ow-k (with a cutoff momentum A = 2 fm^^). It is worth noting that the degree of accuracy 
shown in Table [2] allows to reproduce the /datum of the original potential up to i?iab = 

Table 2: np ^Sq phase shifts (deg) as predicted by the CD-Bonn potential and its 14ow-k (A = 2.0 fm~^). 



E\ah 


CD-Bonn 


^ow-k 


Expt. 


1 


62.1 


62.1 


62.1 


10 


60.0 


60.0 


60.0 


25 


50.9 


50.9 


50.9 


50 


40.5 


40.5 


40.5 


100 


26.4 


26.4 


26.8 


150 


16.3 


16.3 


16.9 


200 


8.3 


8.3 


8.9 


250 


1.6 


1.6 


2.0 


300 


-4.3 


-4.3 


-4.5 



In Table [3], for sake of completeness, we present also the deuteron quadrupole moments for different 
values of the cutoff momentum A, calculated with the corresponding Viow-k's, derived from the Argonne 
Vi8 potential. These values are compared with the one predicted by the original potential. 

Table 3: Deuteron quadrupole moments (fm^) as predicted by the Argonne Vis potential and its Mow-k 
for different values of A (fm~^) 



Argonne Vig 


A 


Mow-k 


0.270 


1.8 


0.268 




2.0 


0.269 




2.2 


0.270 




2.4 


0.270 




2.6 


0.270 




2.8 


0.270 



In Figs. [20] and 1211 we compare np ^Sq T{p\p^p^) matrix elements given by the CD-Bonn potential 
and its Viow-k with a cutoff momentum A = 2.0 fm~^ and p = 0.35, 1.34 fm~^. From these figures it 
can be seen that, even if the half-on-shell T-matrix is not exactly preserved by the ALS transformation, 
it is reproduced to a good degree of accuracy. 

As mentioned earlier, modern A^A^ potentials all possess a strong short-range repulsion. As a result, 
the fc-space matrix elements {k' \ Vnn \ k) are large for large momenta, as illustrated in Fig. [221 where 
we compare the diagonal matrix elements of the phase-shift equivalent potentials Nijmll, Argonne Vig, 
CD-Bonn, and N^LO. 

It is of interest to examine the behavior of Vow-k- In Fig- [231 display some momentum-space 
matrix elements of Vow-k derived from the above V^n^s. Here, we see that Vow-k niatrix elements vary 
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Fig. 20: (Color online) Comparison of the half-on-shell T{p',p,p'^) matrix elements (continuous line) 
with those of Tiov,-k{p' ^P^p"^) for P = 0.35 fm~^ (-Eiab = 10 MeV). The asterisk indicates the on-shell 
matrix element. The cutoff momentum is A = 2 fm~^. 
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Fig. 21: (Color online) As in Fig. [201 but for p = 1.34 fm-^ (Eiab = 150 MeV). 



Table 4: Calculated Pd's with different A^A^ potentials and with the corresponding l^ow-k's- 





Nijmll 


Argonne Vis 


CD-Bonn 


N^LO 


Full potential 


5.63 


5.76 


4.85 


4.51 


Vlow-k 


4.32 


4.37 


4.04 


4.32 



smoothly and become rather small near the cutoff momentum A = 2.1 fm~^. Clearly, Viow~k is a smooth 
potential, no longer having the strong short-range repulsion contained in V^n. 

From Fig. [22] it is evident that the matrix elements of VivAr for the various potentials considered 
are all very different, even if the latter are onshell equivalent. It is interesting to note that the Viow-k's 
derived from them are very close to each other, as shown in Fig. [231 This suggests that the various 
V|ow-k's are almost identical. 




0.0 1.0 2.0 3.0 4.0 

k [fm"'] 



Fig. 22: (Color online) Matrix elements in the ^5*0 channel for different A^A^ potentials. 

In Ref. |163j this V^ow-k property has been evidenced through the study of the ground-state properties 
of some doubly-magic nuclei starting from the above mentioned A^A^ potentials. As a matter of fact, it 
has been found that the calculated binding energies per nucleon and the rms charge radii are scarcely 
sensitive to the choice of the potential from which the Viow-k is derived. This insensitivity may be traced 
to the fact that when renormalizing the short-range repulsion of the various potentials, the differences 
existing between their offshell properties are attenuated. In this regard, let us consider the offshell 
tensor force strength. This is related to the D-state probability of the deuteron Pd, which implies 
that, when comparing A^A^ potentials, offshell differences are seen in Pd differences (see for instance 
Ref. [53]). In Table [H the predicted Pd's for each of the potentials under consideration are reported and 
compared with those calculated with the corresponding l^ow-k's (A = 2.1 fm~^). It can be seen that 
while the Pd's given by the full potentials are substantially different, ranging from 4.5 to 5.8 %, they 
become very similar after renormalization. This is an indication that the starting "onshell equivalent" 
potentials have been made almost "offshell equivalent" by the renormalization procedure. 

It is worth to point out that when the short-range repulsion of different Vattv's is renormalized by 
means of the G-matrix approach the differences between these potentials are strongly quenched too. 
This is illustrated in Fig. [211 There we report the correlation plot between Nijmegen II and CD-Bonn 
matrix elements in the HO basis with the oscillator parameter hu = 14 MeV as well the correlation 
plot between the corresponding Gt matrix elements with a Pauli operator tailored for the sd shell (see 
Sec. 14.1.21) and a starting energy of -5 MeV. We see that while the bare matrix elements are quite 
different this is not the case for the Gt matrix ones. 
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Fig. 23: (Color online) As in Fig. [22l but for the respective Viow-k's with a cutoff momentum A = 2.1 
fm~^. 
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Fig. 24: (Color online) Correlation plots between Nijmegen II and CD-Bonn matrix elements and 
between the corresponding Gt matrix elements (see text for details). The red and blue dots correspond 
to bare Vnn and Gr-niatrix elements, respectively. 
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Before closing this section, a brief discussion of Vfow-k in the framework of the RG-EFT approach 
may be in order. A main point of the RG-EFT approach is that low-energy physics is not sensitive 
to fields beyond a cutoff scale A. Thus for treating low-energy physics, one just integrates out the 
fields beyond A, thereby obtaining a low-energy effective field theory. In RG-EFT, this integration 
out, or decimation, generates an infinite series of counter terms |158] which is a simple power series in 
momentum 

The question as to wether the decimation of high-momentum modes in the derivation of V^ow-k 
generates a series of counter terms has been recently investigated in Refs. [164j and |165] . In these 
studies it is assumed that V^iow-k can be expressed as 

Viow-k(g, q') ^ VNNiq, q') + Kountorlg, q')] (g, q') < , (139) 

where the counter potential is given as a power series 
Kountcr(g, q') = Co + C,q' + C',q" + C,{q' + q") + C',q\" + C,{q' + q'') + C'.q'q" + C'^q'q" + ... (140) 

The counter term coefficients are determined using standard fitting techniques so that the right 
hand side of Eq. (11391) provides a best fit to the left-hand side of the same equation. In Ref. |164j this 
fitting was performed over all partial wave channels and a consistently good agreement was obtained. 

An examination of the counter terms was made in [1641 1165] , where the coefficients of some of them 
are also listed. A main result is that the counter terms are all rather small except for Cq and C2 of the 
S waves. This is consistent with the RG-EFT approach where the counter term potential is given as a 
delta function plus its derivatives [158] . 



5 Realistic shell-model calculations 
5.1 The early period 

As discussed in the Introduction, while the first attempts to derive the shell-model effective interaction 
from the free A^A^ potential go back to some forty years ago, the practical value of this approach has 
emerged only during the last decade or so. On this basis, we may call the "early period" that extending 
over the previous three decades. It is the aim of this section to give a brief survey of the main efforts 
and developments made in the field during this long period of time. 

The first attempt to perform a realistic shell-model calculation may be considered the one by Dawson, 
Talmi and Walecka [17]. They chose ^'^O as testing ground and approximately solved a Bethe-Goldstone 
equation for the wave function of the two valence neutrons using the free A^A^ potential of Brueckner- 
Gammel-Thaler |142] which contains a hard core. This initial approach, in which the effective interaction 
was simply taken to be the Brueckner reaction matrix, was followed by several authors |166l 11671 11681 
1169] . In this context, the merit of the various methods used to obtain the nuclear reaction matrix 
elements was a matter of discussion |170[ I171[ I172^ 1173] . 

A quantitative estimate of the core polarization effects was first given by Bertsch [174] . who found, 
for the Kallio-Kolltveit force [168] , that they can produce a correction to the interaction of two valence 



nucleons in ^^O and ^^Sc as large as 30% of the G-matrix interaction. Soon afterwards a turning point 
in the development of the field was marked by the work of Kuo and Brown who derived an sd- 
shell effective interaction from the Hamada- Johnston potential [63] . After having carefully evaluated 
the reaction matrix elements for this potential, they took into account the corrections arising from 
the one-particle-one-hole {Ip — Ih) excitation of the ^^O core (the so-called "bubble diagram" Gipih). 
Comparison of the calculated and experimental spectra of ^^O and ^^F provided clear evidence of the 
crucial role played by the core polarization. 
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This remarkable achievement prompted the derivation of effective interactions for nuclei in other 
mass regions, such as the calcium [175] , nickel |176i H77\ and lead region [178^ I179t 1180] . These effective 
interactions were all deduced from the Hamada- Johnston potential by use of the G-matrix theory 
and including core-polarization effects. In this context, the individual effects of the various second- 
order diagrams were studied |18H I175j and the conclusion was reached that the inclusion of the core 
polarization diagram, Gipih, is generally sufficient to obtain a satisfactory agreement between calculated 
and experimental spectra |182j . 

The substantial soundness of the effective interactions derived along the lines of the Kuo-Brown 
approach is evidenced by the fact that, with empirical modifications, they have been used in a number 
of successful shell-model calculations to date. This has brought about various "semi-realistic" versions 
of the original Kuo-Brown interaction. In this connection, it is worth pointing out that these empirical 
modifications concern essentially the monopole parts of the effective interaction, as described in detail 
in Ref. |183] . where the KBl, KB2, and KB3 interactions were introduced and applied to the study of 
/p-shell nuclei. This kind of approach, which is still favored by several authors, especially in the context 
of large scale shell model calculations, is outside the scope of the present article. We only mention here 
that recently three-body monopole corrections to realistic forces have been considered as a simple way 
to avoid a full treatment of three-body forces p[84] . For a detailed discussion of the various effective 
interactions and the results of several applications, the reader is referred to the review papers by Brown 
Otsuka et al. [IS], Talmi [S], and Caurier et al. [57], where a full list of references can also be 
found. 

It is interesting to note that in spite of the very encouraging results produced by the first generation 
of realistic shell-model calculations, the general attitude of the workers in the field was rather pessimistic. 
This is well reflected in the final statements made by Kirson in his Summary talk |185] at the 1975 
Tucson International Conference: 

"To summarize my summary, I would say that there are definitely major obstacles in the 
way of doing really convincing calculations. We have a theory, we are aware of the weakness 
of our computations, and we should certainly invest the effort needed to plug the more 
obvious holes. But we should recognize the lack of precision inherent in our inability to 
go to high orders in perturbation theory and adopt more qualitative methods of extracting 
information and gaining understanding. New ideas, new viewpoints are much needed. We 
should also try to withstand the temptations of seductive ^^O, and pay some attention to 
other systems." 

The aim of the Tucson Conference [186j was to define the current status of effective interaction and 
operator theory in the mid 1975s. The main issues of this Conference were expounded in the subsequent 
excellent review by Ellis and Osnes |126] . As it emerges from the first three talks |187tll88|ll89j presented 
at the Tucson Conference and as mentioned in Sec. [3l in the years around 1970 the theoretical bases 
for an energy-independent derivation of the effective interaction were laid down within the framework 
of the time-independent |117j as well as time-dependent [1201 1121] perturbation theory. At the same 
time, the problems inherent in the practical application of the theory became a matter of discussion. 
In this context, considerable effort was made to study the convergence properties of the perturbative 
approach to the derivation of the effective interaction. As representative examples for these studies, 
we may mention here the works of Refs. |190i I19H I192[ I193j . In the first of these works all significant 
terms through third order and a few selected fourth-order terms in the G matrix were calculated and 
the convergence order-by-order of the perturbation series was investigated, while in the second one the 
focus was on the convergence rate of the sum over intermediate-particle states in the second-order core 
polarization contribution to the effective shell-model interaction. In the third work the authors raised 
the problem of the divergence of the perturbation series due to the occurrence of the so-called "intruder 
states" . 
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A most important turn in the field of realistic shell-model calculations resulted from the explicit 
inclusion |194l [T^J of the folded diagrams in the derivation of the effective interaction within the 
framework of the Kuo, Lee, and Ratcliff theory [121] . In these works, starting from a realistic G matrix, 
effective interactions for scZ-shell nuclei were derived by means of the folded-diagram expansion, and 
employed in calculations ioi A = 18 to 24. Note that in Ref. |123j the very convenient method proposed 
in |130] (see Sec. 13. 3p was used to sum up the folded diagram series. It turned out that the contribution 
of the folded diagrams to the effective interaction is quite large, typically about 30%. In the earlier 
study [194] the G matrix was obtained starting from the Reid soft core potential while in the later 
study |123] both this potential and the Paris potential were used. It may be worth noting that this 
latter paper was the first one where the Paris potential was employed in a microscopic calculation of 
finite open-shell nuclei. 

As pointed out in Sec. 12. 1^ the advent of advanced meson-theoretic A^A^ potentials, such as the Paris 
and Bonn ones, was instrumental in reviving the interest in the realistic effective interactions during 
the 1980s. In this context, a main issue was the sensitivity of nuclear structure calculations to the 
strength of the tensor-force component contained in the A^A^ potential, for which a practical measure 
is the predicted D-state probability of the deuteron Po (see Sec. 12.11) . 

Earlier work on this problem was done by Sommermann et al [195J who reexamined the intermediate 
state convergence problem for the core polarization diagram in the expansion of an effective shell-model 
interaction using the Bonn-Jiilich (MDFPA2) potential |196l I197j in comparison to the Reid soft-core 
interaction [M]. It is a remarkable achievement of this work to have shown that the so-called Vary- 
Sauer-Wong effect, namely the slow rate of convergence found for the Reid potential |19H 1198] . was of 
negligible importance for the Bonn-Jiilich potential, as a consequence of the much weaker tensor force 
contained in the latter. 

Some ten years later, using a G-matrix folded-diagram method, the sci-shell effective interaction was 
derived from various potentials and apphed to calculate the spectra of some light S(i-shell nuclei |199j . 
In particular, the three variants of the OBE parametrization of the Bonn full model, Bonn A, Bonn B, 
and Bonn C, differing in the strength of the tensor force (see Sec. 12. ip as well as the Paris potential, 
were considered. It turned out that the best agreement with experiment was achieved with the weakest 
tensor force potential (Bonn A). 

The development of more realistic A^A^ potentials, such as the Bonn and Paris potentials, led also in 
the early 1990s to new studies of the role of the third- and higher-order perturbative contributions to 
the effective interaction, which had been investigated some twenty years before in the work of Ref. [190] . 
Within the framework of the folded diagram theory, the inclusion of the third-order diagrams in G in 
the calculation of the Q-box was studied and the order-by-order convergence of the effective interaction 
in terms of the mass number A was examined |200[ 12011 12021 1203] . The main findings of the last study of 
this series |203] were that the convergence of the perturbation expansion seems to be rather insensitive 
to the mass number and that the effects of third-order contributions in the T = 1 channel are almost 
negligible. We refer to Ref. [45j for a review of the status of the theory of realistic effective interactions 
in the mid 1990s. 

We have given above a brief review of the great deal of work done in the field of realistic shell-model 
calculations during what we have called, somewhat arbitrarily, the "early period" . Before finishing this 
subsection, a word is in order concerning this choice of presentation, which, in our opinion, reflects 
how hard it has been to make this kind of calculations a reliable tool for quantitative nuclear structure 
studies. 

We draw the reader's attention to the fact that, after the calculations performed by Kuo and 
collaborators from the mid 1960's up to the mid 1970's [HSl [HSl [IHl [HSl [173 [ISU!, practically no 
other calculations were done for about two decades having as only aim a detailed description of nuclear 
structure properties. This lack of confidence in the predictive power of realistic effective interactions is 
well evidenced by the fact that in almost all the papers on this subject published from the mid 1970's up 
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to about the early 1990's, the focus has been on nuclei of the srf-shell, considered as a necessary testing 
ground for the effective interaction before embarking on a more ambitious program of applications. In 
this connection, the warning by Kirson (see above citation) against the temptations of seductive ^^O 
acquires a particular significance. 

5.2 Modern calculations 

As discussed in the previous subsection, during the 1980s a substantial progress was made in the 
practical application of shell model with realistic effective interactions. This has resulted in a new 
generation of realistic shell-model calculations oriented to the study of the spectroscopic properties of 
nuclei in various mass regions. 

The first attempts in this direction focused on the Sn isotopes [204^ 133 I2U5] , in particular on the 
neutron deficient ones. In the work of Ref. [204j the Bonn A potential was used while in that of Ref. [15] 
results obtained with the stronger tensor-force versions of this potential (Bonn B and Bonn C) were 
also reported. In both these works all two-body diagrams through third order in the G matrix were 
included in the Q-box. In the study of Andreozzi et al. |205j both the Paris and Bonn A potentials 
were employed and the Q-box was taken to be composed of diagrams through second order in G (a 
comparative discussion of the above studies can be found in f205j). The results of these initial "modern" 
calculations, in particular those obtained with the Bonn A potential, turned out to be in remarkably 
good agreement with the experimental data. 

These achievements encouraged further work along these lines. In the few following years realistic 
shell-model calculations were carried out for a number of nuclei around doubly magic ^°°Sn |134] . ^^^Sn 
[2061 [2071 [2081 [209] . and ^ospb (2l0l M^- Note that in all these calculations, except the one of 
Ref. [209] where the CD-Bonn 1996 potential [87| was employed, the effective interaction was derived 
from the Bonn A potential. We refer to Covello et al. [39J for a survey of the status of realistic shell- 
model calculations at the end of the 1990s. Suffice it to say here that the rms deviation for the spectra 
of 15 medium- and heavy-mass nuclei reported in this paper is in 7 cases less than 100 keV, only in one 
case reaching 165 keV. 

The success achieved by these calculations gave a clear-cut answer to the crucial, long-standing 
problem of how accurate a description of nuclear structure properties can be provided by realistic 
effective interactions, opening finally the way to a more fundamental approach to the nuclear shell 
model than the traditional, empirical one. 

As is well known, much attention is currently being focused on nuclei in the regions of shell closures 
off stability, as they allow to explore for possible shell-structure modifications when approaching the 
proton and neutron drip lines. The experimental study of these nuclei is very difficult, but in recent 
years new data have become available for some of them, which provide a challenging testing ground 
for realistic effective interactions. Starting in the early 2000s, several "exotic" nuclei around doubly 
magic ^''^Sn and ^^^Sn have been studied employing realistic effective interactions derived from the high- 
precision CD-Bonn A^A^ potential [88]. In the first paper of this series |213j shell-model calculations 
were performed for Sn isotopes beyond A^ = 82, in particular the heaviest Sn isotope known to date, 
^^^Sn. In this work, similarly to what had been done in their previous papers, the authors employed a 
G-matrix folded-diagram formalism. A very good agreement was found between the calculated results 
and the available experimental data. 

A most important turn in the derivation of the effective interaction from the free A^A^ potential has 
resulted from the development of the Vjow-k approach (see Sec. 14.21) . This approach to renormalize the 
short-range repulsion of Vat at has been recently used in several calculations. Attention has been focused 
on exotic nuclei neighboring i°°Sn [213], and ^^^Sn [2151 [HSl [HH dSSl [HHl [HH] . Actually, most of these 
works are concerned with nuclei in the latter region, for which there are new interesting experimental 
data. The agreement between theory and experiment has turned out to be very good in all cases. This 
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supports confidence in the predictions of tliese realistic calculations, which may therefore stimulate and 
be helpful to, future experiments. To illustrate this important point, we shall review in Sec. 15.3.31 some 
selected results for ^^^Sn neighbors with neutrons beyond the 82 shell. 

While most of the realistic shell-model calculations carried out in last few years have employed 
effective interactions derived from the CD-Bonn A^A^ potential, some attempts to put to the test chiral 
potentials have also been made. Representative examples for these studies are the works of Refs. [220^ 
12211 1222^ [223| . all of them making use of the ^ow-k renormalization method. Actually, only in |22Uj . 
where the NNLO Idaho potential |lU4j was used, shell model calculations for various two valence-particle 
nuclei were performed. In the other three works [22H 12221 1223j ground-state properties of doubly magic 
nuclei ^He, ^^O, and "^"Ca were calculated within the framework of the Goldstone expansion |122j with 
the N^LO potential [102] . The results obtained in [220] turned out to be comparable to those produced 
by effective interactions derived from the CD-Bonn potential. Clearly, this is related to the important 
question of how much nuclear structure results depend on the NN potential one starts with. We shall 
discuss this point in Sec. I5.3.2[ 

5.3 Results 

5.3.1 Comparison between G-matrix and V^ow-k approaches 

As evidenced by the discussion given in Sees. 15.11 and 15.21 for about forty years the Brueckner G matrix 
represented practically the only tool to renormalize the repulsive core contained in the NN potential. 
After its introduction, the V|ow-k potential has been readily seen as an alternative to the G matrix. In 
fact, V^ow-k is a smooth free-space potential and its matrix elements between unperturbed two-particle 
states can be directly used as input to the Q-box folded-diagram procedure of Sec. [31 Therefore, a main 
issue has been to assess the merit of the Vio„_k approach in practical applications. 

First work along this line was done by Kuo et al. j224j . who carried out shell-model calculations 
for ^^O using two-body matrix elements of effective interactions derived from the Bonn A and Paris 
potentials. Other studies ^225| H5| I226| were then performed which, starting from the CD-Bonn poten- 
tial, extended the comparison between the G matrix and V^ow-k approaches to nuclei heavier than ^'^O. 
In particular, the calculations of Ref. [225j concern ^^^Sn neighbors, while results for ^^O, ^^^Te, and 
^^^I are presented in [33]. A comparison between the G matrix and \^ow-k spectra for the heavy- mass 
nucleus ^^°Po can be found in [226j . In all these calculations the cutoff parameter A is chosen in the 
vicinity of 2 fm~^. 

It has been a remarkable finding of all these studies that the V^ow-k results are as good or even slightly 
better than the G-matrix ones. The rather close similarity between the results of the two approaches 
finds an exception for the Paris potential [224] . We shall comment on this point in Sec. 15.3.21 It is also 
important to point out that the l^ow-k approach appears to work equally well in different mass regions. 
A detailed description of the above mentioned calculations can be found in the cited references. Here, 
for the sake of completeness, we report two examples concerning ^^O and ^^^Te, which are both taken 
from [l3], where more details about the calculations can be found. 

The experimental and calculated spectra of ^^O including the five lowest-lying states, and those of 
^^^Te with the eight lowest-lying states are reported in Figs. [25] and [26] respectively. The calculated 
spectra have been obtained from the CD-Bonn potential through both the Viow-k and G-matrix folded- 
diagram approaches. As for the Mow-k calculations, two values of the cutoff momentum have been used, 
A = 2.0 and 2.2 fm~^. For ^^O we have adopted a model space with the two valence neutrons in the 
(0(^5/2, 0(i3/2, lsi/2) shell, while for ^^^Te the two valence protons are allowed to occupy the levels of 
the (05(7/2, 1(^5/2, l'^3/2, 2si/2, 0/111/2) shell. In both cases the SP energies are taken from experiment, in 
particular from the spectra of ^^O and ^^^Sb for the sd and sdgh shells, respectively. The energy of 
the 2si/2 proton level, which is still missing, is from the study by Andreozzi et al. |206j . Note that in 
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Figs. [25] and [26] the energies are relative to the ground state of the neighboring doubly closed nucleus ^^O 
and ^^^Sn, respectively, with the mass-excess values for ^"^O and ^^^Sb needed for the absolute scaling of 
the SP energies coming from Refs. [227] and [228] . As for ^^^Te, the Coulomb contribution was taken 
into account by adding to the CD-Bonn effective interaction the matrix element of the Coulomb force 
between two gi protons with J = 0. 
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Fig. 25: Low-lying states in ^^O. 

As pointed out above, we see from Figs. [25] and [26] that the l^ow-k and G-matrix results are rather 
close to each other and in quite good agreement with experiment. It is worth noting that the Mow-k 
results for the two values of A do not differ significantly. We shall come back to this point later. 

In this context, we may also mention the work of Refs. [229j and |230j . where harmonic-oscillator 
matrix elements of the Viow-k potential and G matrix are compared in 4 major shells. It was found that 
both the T = and T = 1 matrix elements are very similar. 

On the above grounds, the conclusion is reached that the Mow-k approach is a reliable way to 
renormalize the bare NN potential before employing it in the derivation of the shell-model effective 
interaction. 

Another relevant point to be stressed is that the Viow-k is far easier to be used than the G matrix. 
In fact, Mow-k is a soft NN potential and as such does not depend either on the starting energy or on 
the model space, as is instead the case of the G matrix, which is defined in the nuclear medium. This 
has the desirable consequence that the V^ow-k niatrix element to be used in an interaction vertex of a 
given Q box diagram is completely determined by the incoming and outcoming lines attached to the 
vertex under consideration. Hence, it does not depend on any other part of the diagram. Also, it is 
worth noting that the same l1ow-k can be used in different mass regions. 

Finally, when considering the features of the V|ow-k approach, one has also to keep in mind that, as 
discussed in Sec. 14.2.31 different NN potentials lead to low-momentum potentials which are remarkably 
similar to each other. In the same section, we have showed that the differences between the different 
potentials are also quenched within the G matrix approach. This will be further discussed in Sec. 15.3.21 

We conclude by mentioning that the V^ow-k approach has also proved to be an advantageous alter- 
native to the G matrix in the calculation of the ground-state properties of doubly closed nuclei. As 
is well known, the traditional approach to this problem is the Brueckner-Hartree-Fock (BHF) theory. 
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In the study of Ref. [2'21\ it has been showed that the Viow-k is suitable for being used directly in the 
Goldstone expansion. Good results for ^^O and ^°Ca were obtained by solving the Hartree-Fock (HF) 
equations for V^ow-k and then using the resulting self-consistent basis to compute higher-order Goldstone 
diagrams. In this connection, it should be mentioned also the work of Ref. [231] , where use was made 
of the Viow-k potential in a HF and BHF approach to evaluate the binding energy and the saturation 
density of nuclear matter and finite nuclei. Their conclusion was that for nuclear matter at small and 
medium densities as well as for finite nuclei HF and BHF results with \4ow-k are nearly identical. 
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Fig. 26: Low- lying states in ^^^Te. 



5.3.2 Calculations with different A^A^ potentials 

As discussed in Sec. O a main problem relevant to realistic shell-model calculations has been to assess 
the accuracy of the effective interaction derived from the A^A^ potential. However, it is also an important 
issue what information comes out on the A^A^ potential from this kind of calculations. As we have seen 
in Sec. [21 from the early 1950s on several different potentials have been constructed. The potentials 
of the last generation all fit equally well the A^A^ phase shifts, but, except for the OPE tail, they 
are quite different from each other. Actually, these potentials are well-constrained by the deuteron 
properties and the A^A^ scattering data for i^iab < 350 MeV and therefore their offshell behavior is 
subject to ambiguities. This is reflected in the different values predicted for the D-state probability of 
the deuteron Pd (see Sec. [2]) as well as in the different binding energies calculated for A = 3,4 systems 
[232] . In fact, for A > 2 the interaction may be off the energy shell since two nucleons may have different 
energies before and after they interact. 

In this context, it is of great interest to perform nuclear structure calculations to see to which extent 
the differences between A^A^ potentials persist once the short-range repulsion is renormalized through 
the G-matrix or the l^ow-k approach. This problem has long been investigated by studying properties 
of both nuclear matter and doubly closed nuclei. A list of references on this subject can be found in 
the work of Ref. p22j , which we have briefly discussed in Sees. 14.2.31 and 15. 2[ As for open-shell nuclei. 



some attempts have been made within the framework of the shell model. 
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In the early work by Sommermann et al. [195J , the comparison between Bonn- Jiihch and Reid resuhs 
aimed at investigating the intermediate-state convergence problem, as discussed in Sec. I5.1[ The paper 
by Shurpin et al. [123J had instead as a main purpose to compare the spectra of several sd-shell nuclei 
calculated with two different potentials, Paris and Reid. The results showed that the Paris potential 
was as good as, if not better than, the Reid potential. The idea of comparing shell-model calculations 
with different NN potentials was revived at the end of the 1980s. In the works of Refs. |233[ 1199] and 
[234] sci-shell nuclei were studied using as input several different NN potentials, all constructed before 
the 1990s. These works were based on the use of the reaction matrix G to derive the effective interaction 
through the Q-box folded- diagram method, including up to second-order diagrams. The outcome of 
these studies was that the effective interaction matrix elements become more attractive when a A^A^ 
potential with a weaker tensor-force strength is used. Jiang et al. [199] noted that the nonlocality of 
the potential gives a further contribution in this direction. All these calculations indicated in general 
that the best agreement with experiment was produced by the Bonn potential. 

A confirmation of the important role played by weaker tensor force potentials came out from the 
paper of Ref. [205] . where Paris and Bonn A results for neutron deficient Sn isotopes are compared. A 
better agreement with experiment was obtained with the Bonn A potential. The merit of potentials 
with weak tensor force was also evidenced in [45], who performed realistic shell-model calculations for 
nuclei in the mass regions of ^^O, ^"Ca, and ^°°Sn using the three variants (A, B, and C) of the Bonn 
potential. 

A comparative study of shell-model results obtained from different A^A^ potentials, including the 
high-precision CD-Bonn potential, was made for the first time in |235j . In this paper, the nucleus ^^^Te 
was used as a testing ground to investigate the effects of the CD-Bonn potential with respect to those 
of the three potentials Bonn A, Paris, and Nijm93. As in all previous calculations, it was found that a 
weak tensor force leads to a better description of nuclear structure properties. Furthermore, it turned 
out that the calculated spectra, except that obtained from the Paris potential, are only slightly different 
from each other. 

All the above calculations made use of the reaction matrix G and their results provide evidence 
that the differences existing between various potentials are quenched by the G-matrix renormalization 
procedure (see Fig. [21]). One may wonder what happens when the V^ow-k potential is used. The first 
answer to this question may be found in |224j . where the Viow-k approach is used to calculate the Paris 
and Bonn A spectra for ^^O. It is interesting to note that these two spectra are much closer to each 
other than those obtained with the G matrix approach |199] . This reflects the fact that the Viow-k 
procedure is more effective in reducing the differences between the two potentials than the G-matrix 
approach. The matrix elements in the ^Si channel for both Paris and Bonn A potentials are reported 
in Fig. 1 of [224] , where it is shown that they become similar after the Viow-k renormalization, the latter 
being much larger for the strongly repulsive Paris potential. 

In Sec. 14.2.31 (see Figs. [221 and [23]) . we have reported a similar comparison in the ^5*0 channel for 
the chiral potential N^LO and the three high-precision potentials CD-Bonn, Nijmll, and Argonne Vi^. 
The differences between the matrix elements of these potentials practically disappear when the V|ow-k 
matrix elements are considered. On these grounds one expects that shell-model effective interactions 
derived from phase-shift equivalent A^A^ potentials through the Viow-k approach should lead to very 
similar results. 

This we have verified by calculating the spectra of ^^"^Te starting from the CD-Bonn, Nijmll, and 
Argonne Vis potentials. For all the three V^ow-k's derived from these A^A^ potentials, use has been 
made of the same cutoff momentum A = 2.2 fm~^. Other details on how these calculations have been 
performed may be found in Sec. 15.3.31 where we present results for nuclei beyond ^^^Sn, which have 
been obtained using the CD-Bonn potential. 

In Fig. [271 fhe three calculated spectra are shown and compared with the experimental one. Note 
that, with respect to Fig. [26l the energies are relative to the ground state and higher lying levels are 
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Fig. 27: Spectrum of ^^^Te. Predictions by various A^A^ potentials are compared with experiment. 
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Fig. 28: Experimental and calculated spectra of ^^^Sn. 



included. We see that the calculated spectra are very similar, the differences between the predicted 
energies not exceeding 80 keV. It is also seen that the agreement with experiment is very good for all 
the three potentials. In fact, the rms deviation is 115, 128, and 143 keV for CD-Bonn, Argonne Vig, 
and Nijmll, respectively. 

Finally, it is of interest to make a few comments regarding the value of the cutoff momentum A. 
In Sec. 14.2.31 the plausibility of using a value in the vicinity of 2 fm^^ was discussed. It was also 
shown in Sec. 15.3.11 that shell-model results do not change significantly when varying A from 2 to 2.2 
fm~^. A preliminary study of the dependence of shell-model results on A can be found in the work of 
Ref. [236] ■ where results for ^^^Te are reported using the CD-Bonn, Argonne Vis, and Nijmll potentials 
with A varying over a rather large range (1.5-2.5 fm~^). It has been found that the best agreement 
with experiment for each of these three potentials corresponds to a somewhat different value of A. 
However, changes around this value do not significantly modify the quality of agreement for all the 
three potentials. As an example, the rms deviation for Nijmll remains below 150 keV for any value of 
A between 1.8 and 2.2 fm~^ This subject deserves certainly further studies. 

5.3.3 Neutron rich nuclei beyond ^^^Sn: comparing theory and experiment 

The study of exotic nuclei is a subject of a current experimental and theoretical interest. Experimental 
information for nuclei in the vicinity of ^^Ni, ^°°Sn, ^'^^Sn, which have been long inaccessible to spec- 
troscopic studies, is now available thanks to new advanced facilities and techniques. The advent of 
radioactive ion beams is opening new perspectives in this kind of studies. Nuclei in the regions of shell 
closures are a source of direct information on shell structure. New data offer therefore the opportunity 
to test the shell model and look for a possible evolution of shell structure when going toward proton or 
neutron drip lines. 

This is stimulating shell-model studies in these regions. Several realistic shell-model calculations 
have been recently performed for nuclei around ^''^Sn and ^^^Sn, as mentioned in Sec. 15. 2[ We report 
here some selected results on ^^^Sn neighbors with > 82. More precisely, we consider the three nuclei 
^^^Sn, ^^^Sb, and ^^^Sb, which have been very recently studied in [219] . [218j . and |163] . to which we 
refer for a more detailed discusssion. 

The experimental data which have been now become available for these neutron-rich nuclei may 
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Fig. 29: Experimental and calculated spectra of ^^^Sb. 



suggest a modification in the shell structure. They are, in fact, somewhat different from what one 
might expect by extrapolating the existing results for < 82, and as a possible explanation a change 
in the single-proton level scheme has been suggested. The latter, caused by a more diffuse nuclear 
surface, may be seen as a precursor of major effects which should show up at larger neutron excess. 

We shall see now that the properties of these three nuclei are all well accounted for by a shell- 
model calculation with a two-body effective interaction derived from a modern NN potential. In so 
doing, we evidence the merit of realistic shell-model calculations, showing, at the same time, that 
there is no need to invoke new effects to describe few-valence-particle nuclei just above ^^^Sn. In this 
connection, it is worth mentioning that different shell-model calculations have been performed for the 
above nuclei |237l I238[ I239[ 12401 12411 1242] . None of these studies, however, has succeeded in accounting 
simultaneously for the peculiar features of all the three nuclei. 

In our calculations, we assume that ^^^Sn is a closed core and let the valence neutrons occupy the six 
levels Ohg/2, I/7/2, I/5/2! 2^3/2, 2pi/2, and 0^13/2 of the 82-126 shell, while for the proton the model space 
includes the five levels 0*77/2, 1(^5/2, 1*^3/2, 2si/2, and 0/111/2 of the 50-82 shell. The proton and neutron 
single-particle energies have been taken from the experimental spectra of ^^^Sb and ^^'^Sn, respectively, 
the only exceptions being those of the Si/2 proton and ^13/2 neutron levels which are still missing. The 
energy of the former, as mentioned in Sec. 15.3. H is from [206] . while that of the latter from |213j . 
The two-body matrix elements of the effective interaction have been derived from the CD-Bonn A^A^ 
potential. This potential is renormalized by constructing the ^ow-k with a cutoff momentum A =2.2 
fm~^. The effective interaction has been then derived within the framework of the Q-box plus folded 
diagram method including diagrams up to second order in V^ow-k- The computation of these diagrams 
is performed within the harmonic-oscillator basis using intermediate states composed of all possible 
hole states and particle states restricted to the five shells above the Fermi surface. This guarantees 
the stability of the results when increasing the number of intermediate particle states. The oscillator 
parameter used is hu = 7.88 MeV and for protons the Coulomb force has been explicitly added to the 
Viow-k potential. 

The experimental and calculated spectra of ^'^^Sn, ^'^^Sb, and ^'^^Sb are compared in Figs. EHl [291 
and [30], respectively. Note that both the calculated and experimental spectra of ^^^Sb include all levels 
lying below 1 MeV excitation energy, while for ^^^Sb all the yrast states have been reported. As for 
^^^Sn, only the levels reported in Fig. [28] have been observed, while in the calculated spectrum five 
states, predicted in between the 6+ and 8"*" states, are omitted. 

The two isotones with A = 134 are the most appropriate systems to study the effects of the effective 
interaction. In fact, ^^^Sn is a direct source of information on the neutron-neutron channel while ^^^Sb 
on the proton-neutron one. The nucleus ^^^Sb represents a further interesting test, since it allows to 
investigate at once the role of the effective interaction in both channels. 

From the above figures, we see that the experimental levels in all the three nuclei are very well 
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Fig. 30: Experimental and calculated spectra of ^'^^Sb. 



reproduced by the theory, the discrepancies being of the order of few tens of keV for most of the states. 
Note that the very low-energy positions of both the first-excited 2+ and 1~ states in ^^^Sn and ^^^Sb, 
respectively, are well accounted for, as is also the case for the low excitation energy of the 5/2"*" state in 
^^^Sb. We refer to the above mentioned papers for more details on each of the three calculations and for 
a discussion of the structure of the states in ^'^^Sb and ^^^Sb. As for the latter, it is shown in [163] that 
it is the admixed nature of the 5/2+ state that explains its anomalously low position and the strongly 
hindered Ml transition to the ground state. The role of the effective interaction as well as of the 5/2+ 
single-proton energy in determining this situation has been examined, with particular attention to the 
effects induced by the proton-neutron interaction. The relevant role of the latter was also evidenced in 
^^^Sb |218j . For this nucleus, a direct relation can be established between the diagonal matrix elements 
of the proton-neutron effective interaction and the calculated energies reported in Fig. [291 since the 
corresponding wave functions are characterized by very little configuration mixing. In particular, the 
eight lowest-lying states can be identified with the eight members of the 7r(yf7/2 1^/7/2 multiplet while 
the two higher states, with J'" = 1~ and 2~, are the only members of T^d^/2^f7/2 multiplet having an 
experimental counterpart. 

In this connection, it may be interesting to summarize here the analysis of the proton-neutron matrix 
elements reported in the work of Ref. |218j . Focusing on the T^g7/2^f7/2 configuration, we examine the 
various terms which contribute to the effective interaction, namely the Viow-k and all the second-order 
two-body diagrams composing the Q box. It is worth pointing out that the effective interaction 
is obtained by summing to the Q-hox the folded-diagram series [see Eq. (|58|) ]. but the contribution 
of the latter provides a common attenuation for all matrix elements, thus not affecting the following 
discussion. 

The diagonal matrix elements of V|ow-k for the T^gj/2^f7/2 configuration are reported in Fig. ED as a 
function of J together with the contributions arising from the Ip — Ih, 2p — 2h, and ladder excitations. 
In Fig. [32] we present for the same configuration the diagonal matrix elements of the effective interaction 
which we have obtained starting from this l^ow-k- Note that the behavior of V^g shown in Fig. [32]is quite 
similar to that of the T^gT/2^f7/2 multiplet in ^^^Sb, whose members, as mentioned above, are almost 
pure. This behavior, however, is quite different from that of l^ow-k, especially as regards the 0~ — 1^ 
spacing. This means that the second-order two-body contributions accounting for configurations left 
out of the chosen model space are of key importance to renormalize the Viow-k and get the correct 
behavior. In particular, a crucial role is played by the core-polarization effects arising from the Ip — Ih 
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Fig. 31: Diagonal matrix elements of the Viow-k and contributions from the two-body second-order 
diagrams for the 71(77/2^^/7/2 configuration. 

excitations. 

At this point one may wonder why our effective interaction is able to explain the properties of these 
neutron-rich nuclei above ^^^Sn while the realistic shell-model calculations of Refs. |237[ I238[ I239j do 
not. In the latter calculations the effective interaction has been derived within the G-matrix folded- 
diagram approach including in the Q-box diagrams up to third order in G and intermediate states with 
at most 2huj excitation energy [203] . With this interaction, a satisfactory agreement with experiment 
was obtained only for ^^^Sn when using experimental SP energies. A downshift of the vrc?5/2 level with 
respect to the one by 300 keV turned out to be necessary to obtain a good description of ^^^Sb 

[237[ 1239] ■ but did not help for ^^^Sb [238j . We have verified that the differences between the results 
presented in this section and those of Refs. [237112381 1239] can be traced mainly to the different dimension 
of the intermediate state space. In fact, we have found that including intermediate states only up to 
2hw excitation energy in our derivation of the effective interaction, the results of the two calculations 
become very similar. 

In summary, we have shown that modern realistic shell-model calculations are able to describe with 
quantitative accuracy the spectroscopic properties of ^^^Sn, ^^^Sb, and ^^^Sb, in particular, their peculiar 
features. This leads to the conclusion that a consistent description of ^^^Sn neighbors with neutrons 
beyond the N = 82 shell is possible. In this connection, it is worth mentioning two more recent works 
[243[ 1244] . where shell-model results for other nuclei around doubly magic ^^^Sn are reported, including 
predictions for ^^^Sn. 

5.3.4 The role of many-body forces 

In this subsection, we shall discuss the role that the 3A^ forces, both real and effective, play in shell- 
model calculations based on the results we obtain for the even-mass = 82 isotonic chain. This study 
has been performed with the same input for all nuclei of the isotonic chain, namely the SP energies 
and the matrix elements of the two-body effective interaction employed in the calculations discussed in 
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Fig. 32: Diagonal matrix elements of the two-body effective interaction for the Tig-j/2^f7/2 configuration. 



Sec. 15.3.31 Here, the shell-model calculations have been carried out by using the Oslo code 
Let us first consider the two valence-proton system ^^"^Te. 



As mentioned in Sec. 12.31 a genuine 3iV 
potential, which may affect the one- and two-body components of the shell-model effective interaction 
of Eq. (1S5]) . has never been considered explicitly. The contributions of the 3A^ potential to the one- 
body term of H^^ may be assumed as taken into account when, as usual, the SP energies are taken from 
the experimental data. As regards the contribution of the 3iV potential to the two-body component of 
i^gfj, it may be reduced with an appropriate choice of the cutoff momentum A of the l^ow-k- In |229j . 
it was suggested to choose A around 2 fm~^, since this value seems to minimize the role of the 3A^ 
potential when calculating the and ^He binding energies. In this context, it is worth recalling that 
here we have used experimental SP energies and an effective interaction derived from a V|ow-k with 
A = 2.2 fm-^ 

For the ground-state energy of ^^^Te, relative to the doubly magic ^^^Sn, we obtain a value of -20.796 
MeV to be compared with the experimental one -20.584 MeV. As regards the energy spectrum, the rms 
deviation is 0.115 MeV (see Fig. [27j). These results seem to validate for many-nucleons systems the 
hypothesis made in [229j about the choice of A. 
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Let us now come to nuclei with more than two valence nucleons. We know that in this case many- 
body terms appear in the shell-model effective interaction, even if the original potential is only a 
two-body one. The effects of these terms have not yet been quantitatively evaluated pTl [21^ I247j . 
However, it is generally believed that effective many-body forces have little influence on the excitation 
energies, while the binding energies may be significantly affected pL2Y j. In this connection, it may be 
mentioned that Talmi [248] has shown that the experimental binding energies of the tin isotopes can 
be reproduced quite accurately employing an empirical two- and three-body interaction. 

In order to study the effective many-body forces, we have calculated the excitation energies of the 
lowest-lying yrast states in the even = 82 isotones. In Figs. [221 [211 and [23 the calculated energy 
of the J'^ = 2+, 4"'', and 6"^ states, respectively, are reported as a function of the mass number A and 
compared with the experimental data. It is seen that the agreement between theory and experiment is 
of the same quality all over the isotonic chain, thus confirming that effective many-body forces play a 
negligible role in the description of the relative energy spectra. 
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Fig. 34: Same as in Figj33]but for J = 4+. 

In Fig. [36l we report the ground-state energy per valence nucleon. It appears that the experimental 
and theoretical behaviors diverge when increasing the number of protons. As in the case of tin isotopes 
|249j , the binding energy is increasingly overestimated. This deficiency may be healed by adding a small 
repulsive monopole contribution, whose origin may be traced to effective many-body forces [248j. 



6 Summary and conclusions 

In this paper, we have tried to give a self-contained review of the present status of shell-model calcula- 
tions employing two-body effective interactions derived from the free nucleon-nucleon potential. A main 
feature of this kind of calculations, which are commonly referred to as realistic shell-model calculations, 
is that no adjustable parameter appears in the determination of the effective interaction. This removes 
the uncertainty inherent in the traditional use of empirical interactions. 

As discussed in the Introduction, only in the last few years has this approach, which started in 
the 1960s, acquired a reputation for being able to provide an accurate description of nuclear structure 
properties. This means that to place in its proper perspective the current status of the field it is helpful 
to look back and survey the difficulties faced and the progress achieved over about four decades. The 
present review, while aimed at giving an insight into modern realistic shell-model calculations, is also 
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Fig. 36: Ground-state energies per valence proton of = 82 isotones as a function of the valence-proton 
number Z . 
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concerned with the above aspects in regard to both the A^A^ potential and the many-body methods for 
deriving the effective interaction. 

As regards the first issue, in Sec. [2] we have discussed in some detail the main advances in the 
understanding of the A^A^ interaction from the early stages of development until the current efforts 
based on chiral effective field theory. This we have done to show how the advances in this field have 
been intimately related to the development of realistic nuclear structure calculations. The present stage 
of development is certainly enough to provide a sound starting point for a quantitative description of 
nuclear structure properties. It remains a main task for the near future to assess the role of the three- 
body forces in many-nucleon systems. A brief discussion on this matter is given in Sec. 15.3.41 based on 
the results for the even A^ = 82 isotones. 

In Sec. E] we have described how the shell-model effective interaction can be derived within the 
framework of degenerate time- dependent perturbation theory. 

The problem of the short-range repulsion contained in the A^A^ potential has been discussed in 
Sec. m where, after an outline of the traditional Brueckner G-matrix method, we have described a new 
approach to the renormalization of Vnn, which consists in constructing a low- momentum A^A^ potential, 
V|ow-k, that preserves the physics of V^n up to a cutoff momentum A. Since Viow-k is already a smooth 
potential, it can be used directly as input for the derivation of the shell-model effective interaction 
instead of the usual G matrix. While its merit in this respect has been assessed by several calculations 
(Sec. l5.3.Tl) the A^A^ potential l^ow-k is currently attracting much attention. On the one hand, it is being 
applied in various contexts, for instance in the study of few-body problems |250l I114j . On the other 
hand, its properties are being studied in detail, in particular as regards the dependence on the original 
potential model |25H 1222^ 1236] . Concerning shell- model calculations, results obtained with different 
A^A^ potentials using the V^ow-k approach have been presented in Sec. I5.3.2[ 

A main motivation for this review has been the practical success achieved by realistic shell-model 
calculations in the last few years, after a long period of hope and dismay. In Sec. O we have first 
surveyed the main efforts in this field going from the early to present days. We have then presented, by 
way of illustration, some selected results of recent calculations for nuclei around doubly magic ^^^Sn. 
These neutron-rich nuclei, which lie well away from the valley of stability, offer the opportunity for a 
stringent test of the matrix elements of the effective interaction. The very good agreement with the 
available experimental data shown in Sec. 15.3.31 supports confidence in the predictive power of realistic 
shell-model calculations in the regions of shell closures off stability, which are of great current interest. 

As already pointed out above, in this review our focus has been on shell-model calculations employing 
genuine free A^A^ potentials. As a consequence of this choice, even a brief summary of the many successful 
shell-model calculations employing modified versions of free A^A^ potentials would have been beyond 
the scope of the present article. Here, we have also not touched upon another important aspect of shell- 
model work, namely the development of new methods, such as the Monte Carlo shell model (MCSM), 
as well as of high-quality codes, such as ANTOINE |252j . which has greatly extended the feasibility 
of large-scale calculations. However, both these issues are discussed in detail in the recent article by 
Caurier et al. [17] while for a specialized review of the MCSM we refer to Ref. [18] . 

To conclude, we may say that the stage for realistic shell- model calculations is by now well set, 
which is a major step towards the understanding of the properties of complex nuclei in terms of the 
forces among nucleons. In this context, the Viow-k approach provides a new effective tool to handle the 
short-range repulsion of the free A^A^ potential. 
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